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FOREWORD 


This  report  documents  work  accomplished  during  the  period  18 
January  1971  through  17  January  1972  by  Booz,  Allen  Applied  Research, 
Inc.,  P.O.  Box  1797,  Eglin  Air  Force  Base,  Florida,  under  Contract 
F08635-71-C-0093  with  the  U.  S.  Air  Force  Armament  Laboratory,  Eglin 
Air  Force  Base,  Florida.  The  program  monitors  for  the  Armament 
Laboratory  were  Mr.  Dan  Mclnnis  (DLYW)  and  Mr.  Jerry  Bass  (DLYW). 

The  report  consists  of  two  volumes.  Volume  I,  the  User  Manual, 
contains  a  detailed  description  of  the  input  variables,  instructions  for 
placing  the  input  variables  on  punch  cards,  instructions  for  arranging 
the  punch  cards  in  proper  order,  descriptions  and  definitions  of  the  pro¬ 
gram  output,  a  description  of  a  sample  case  which  can  be  used  to  verify 
that  the  program  is  operating  properly,  and  a  checklist  of  input  data 
required  to  make  up  a  test  case.  Volume  II,  the  Analyst  Manual,  con¬ 
tains  the  mathematical  relationships  which  were  used  to  develop  the 
model,  the  assumptions  employed  in  constructing  the  model,  flowcharts 
depicting  the  logical  structure  of  the  model,  a  listing  of  the  source  deck, 
and  a  detailed  description  of  the  simulation  model  based  on  the  coding 
language  employed  and  including  comments  identifying  algebraic  ex¬ 
pressions  which  appear  in  the  mathematical  model. 

The  authors  wish  to  acknowledge  the  assistance  provided  by  several 
individuals  in  developing  the  methodology  and  computer  program  described 
in  this  report.  Mr.  Gary  M.  Grann  of  the  U.  S.  Air  Force  Armament 
Laboratory  was  instrumental  in  developing  the  initial  methodology  and 
computer  program.  Al§o,  the  suggestions  and  comments  provided  by  the 
members  of  the  Methodology  and  Evaluation  Working  Group  of  the 
Degradation  Effects  Program  were  extremely  helpful  in  developing  the 
computer  program. 

This  technical  report  has  been  reviewed  and  is  approved. 


/(IvWUa)  '  i^AAA^/JL 
THOMAS  P.  CHRISTIE 


Chief,  Weapon  Systems  Analysis  Division 


ABSTRACT 


The  Burst  Height  Distribution  (BHD)  Program  described  in  this 
report  was  designed  to  compute  and  display  burst  height  distributions 
for  munitions  aerially  delivered  into  forest  environments.  The  program 
uses  as  input  the  source  and  terminal  X,  Y,  and  Z  coordinates  and  the 
average  diameters  of  branches  surveyed  at  actual  forested  sites,  and 
the  munitions  travel  along  straight-line  trajectories  which  are  randomly 
selected.  Burst  heights  are  computed  for  those  trajectories  which  en¬ 
counter  branches  large  enough  to  detonate  the  munition,  and  after  400 
trajectories  are  examined  (100  from  each  of  four  azimuth  angles),  the 
cumulative  burst  height  distribution  for  the  munition  and  elevation  angle 
is  computed,  printed,  and  optionally  punched  as  output.  The  computer 
program  was  specifically  designed  for  the  Control  Data  Corporation  6600 
computer  system  at  Eglin  Air  Force  Base,  Florida. 


Distribution  limited  to  U.  S.  Government  agencies  only; 
this  report  documents  test  and  evaluation;  distribution 
limitation  applied  January  1072  .  Other  requests  for 
this  document  must  be  referred  to  the  Air  Force  Armament 
Laboratory  (DLYW)  ,  Eglin  Air  Force  Base,  Florida  32542. 


iii 


(The  reverse  of  this  page  is  blank) 


TABLE  OK  CONTENTS 


Section  Page 

I.  INTRODUCTION  .  1 

II.  MATHEMATICAL  MODEL .  3 

Determination  of  Initial  Trajectory  Entry 
Coordinates  Into  Top  Layer .  3 

Translation  and  Rotation  of  Branch 

Coordinates  .  5 

Determination  of  Subsequent  Trajectory 

Exit  and  Reentry  Coordinates .  10 

Determination  of  Munition/Branch 

Encounters  and  Burst  Heights .  14 

Case  1 .  14 

C  ase  2 .  17 

Case  3 .  If) 

Computation  of  Cumulative  Burst  Height 
Distributions .  21 

III.  FLOWCHARTS  .  22 

IV.  SOURCE  LISTING  .  3.3 

V.  SIMULATION  MODEL .  44 

Tables  and  Arrays .  44 

Discussion  of  Simulation .  45 

Program  MAIN .  45 

Subroutine  CFRQ .  57 


v 


LIST  OF  FIGURES 


Figure  Title  Page 

1.  Pictorial  Representation  of  Layers  Used  in 

Forest  Environment  and  Trajectory  Paths .  4 

2.  Relationship  of  Original  and  Transformed 
Branch  Coordinates  After  First  Rotation  in 

X-Y  Plane  Through  Azimuth  Angle  0 .  6 

3.  Relationship  of  Transformed  x',  Y',  and  z! 

Coordinates  to  Xsj,  YSj,  and  Zsj  Coordinates 
After  Rotation  of  Axes  m  Xs-Za  Plane 

Through  Elevation  Angle  cp .  8 

4.  Relationship  of  Entry  and  Exit  Coordinates  as 

Trajectories  Pass  Through  Layers  of  Branches .  11 

5.  Relationship  of  Original  and  Transformed 
Munition  Trajectory  Entry  Point  Coordinates 
After  Rotation  in  X-Y  Plane  Through 

Elevation  Angle  cp .  13 

6.  Munition/Branch  Encounter  -  Case  1 .  15 

7.  Munition/Branch  Encounter  -  Case  2 .  18 

8.  Munition/Branch  Encounter  -  Case  3 .  20 

9.  Flowchart,  Program  MAIN .  23 

10.  Flowchart,  Subroutine  CFRQ .  31 

11.  BHD  Program  Source  Listing 


34 


Tabic* 

I. 


LIST  OF  TABLES 


Title 

Translation  Values  Xt  and  Yt  Used  in 
BHD  Program . 


Page 


0 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS  -  MATHEMATICAL  MODEL 


SYMBOL  OR 
ABBREVIATION 

EQUIVALENT 

IN 

SIMULATION 

MODEL 

DEFINITION 

UNITS 

A 

A 

The  distance  between  the 
branch  ends  along  the  trans¬ 
formed  Z;  axis. 

mm 

B 

B 

The  distance  between  the 
branch  ends  along  the  trans¬ 
formed  Y'  axis. 

mm 

C 

C 

A  value  used  in  computing 
the  distance  in  the  trans¬ 
formed  Y'-Z'  plane  between 
the  munition  trajectory  and 
the  branch  centerline  for 

Cases  2  and  3  [  i.  e. , 

C  --  -Z'jjlB)  -  Y'jjfA)]. 

mm 

Ck 

CFREQU50) 

The  cumulative  frequency 
for  the  k**1  interval  of  the 
burst  height  distribution 
table. 

percent 

Da 

D 

The  square  of  the  distance  in 
tiie  transformed  Y'-Z7  plane 
between  the  munition  tra¬ 
jectory  and  the  branch  center¬ 
line  (see  Figures  6,  7,  and  8). 

mm‘ 

Dd 

BBDIAM 

Tiie  minimum  branch  diam¬ 
eter  required  for  detonation 
by  tiie  munition  body. 

mm 

Fk 

FREQ(1 50) 

The  numbe  ’  of  detonations 
in  the  k^  burst  height  in¬ 
terval. 

none 
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(CONTINUED) 


EQUIVALENT 

IN 

SNMUOi.OR  SIMULATION 
ABBREVIATION  MODEL  DEFINITION 


P 


q 


K 


»>J 


R. 


R 


rn 


^max 


R 


min 


RN 


The  sequential  number  of  the 
rectangle  along  the  X  axis 
defining  the  munition  trajec¬ 
tory  entry  point  into  the  top 
layer  of  branches. 

The  sequential  number  of  the 
rectangle  along  the  Y  axis 
defining  the  munition  trajec¬ 
tory  entry  point  into  the  top 
layer  of  branches. 

AVRAD  The  average  radius  of  the 
branch. 

FRAD  The  radius  of  the  munition 

fuze . 

BRAD  The  radius  of  the  munition 

body . 

The  greater  of  the  distances 
in  the  transformed  -7,' 
plant'  between  the  munition 
trajectory  and  each  end  of 
the  branch. 

The  smaller  of  the  distances 
in  the  transformed  Y'-Z' 
plane  between  the  munition 
trajectory  and  each  end  of 
the  branch. 

A  random  number  selected 
from  a  uniform  distribution 
with  0.  0  <  Rjq  <1.0. 


UNITS 

none 


none 


mm 


mm 


mm 


mm 


mm 


none 


ix 


F 


I 


r 


s. 

p 
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(CONTINU  ED) 

EQUIVALENT 

IN 

SYMBOL  OR  SIMULATION 
ABBREVIATION  MODEL  DEFINITION 


R 


Xi 


RYi 


Hi 


Ra 


Si 


A  t  arn  lorn  integer  between  1 
and  10  selected  from  a  uni¬ 
form  distribution. 

A  random  integer  between  1 
and  10  selected  from  a  uni¬ 
form  distribution. 

[R(l)]i  The  distance  in  the  trans¬ 
formed  Y'-Z'  plane  between 
the  munition  trajectory  and 
one  of  l  he  ends  of  the  branch 
(see  !  igures  7  and  8). 

[R(2)]£  The  distance  in  the  trans¬ 
formed  plane  between 

the  munition  trajectory  and 
one  of  the  ends  of  the  branch 
(see  Figures  7  and  8). 

SI  The  distance  in  the  trans¬ 

formed  Y  -z!  plane  along  the 
center  line  of  the  branch  and 
between  one  of  the  ends  of 
the  branch  and  the  point  at 
which  the  normal  from  the 
munition  trajectory  inter¬ 
sects  the  branch  centerline 
(see  Figure  7). 


MODEL  .  » 


UNITS 

none 


none 


mm 


mm 


mm 


x 
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EQUIVALENT  s 

IN 

SYMBOL  OR  SIMULATION 
ABBREVIATION  MODEL  DEFINITION 


S2 


X 


Ei 


X 


.1 


^max 


S2  The  distance  in  the  trans¬ 

formed  X'-z!  plane  along  the 
centerline  of  the  branch  and 
between  one  of  the  ends  of 
the  branch  and  the  point  at 
which  the  normal  from  the 
munition  trajectory  inter¬ 
sects  the  branch  centerline 
(see  Figure  7). 


XDET  The  transformed  X'  coordi¬ 

nate  of  the  detonation  for 
the  ifh  munition  trajectory. 

XTRAJ(IOO)  The  X  coordinate  of  the  ith 
munition  trajectory  entry 
point  into  the  layer  of 
branches. 


The  transformed  X1  coordi¬ 
nate  of  the  i*n  munition 
trajectory  entry  point  into 
the  layer  of  branches. 

BRANCH(6)  The  X  coordinate  of  an  end 
of  the  jth  branch. 

Xl(1500)  The  transformed  X'  coordi- 
X2(1500)  nate  of  an  end  of  the  j**1 
branch . 


XLIMIT  The  X  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment. 


MODEL 


UNITS 

mm 


mm 


mm 


mm 


mm 


mm 


mm 
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SYMBOL  OR 
ABBREVIATION 

EQUIVALENT 

IN 

SIMULATION 

MODEL 

DEFINITION 

UNITS 

x  min 

1  ho  smaller  of  the  trans- 
forrm  '1  X'  coordinates  of 
l  lie  branch  ends . 

mm 

xOi 

XOUT 

The  X  coordinate  of  the  i^*1 
munit  on  trajectory  exit 
point  from  the  layer  of 
branches. 

mm 

Xri 

The  X  coordinate  of  the  i^ 
munition  trajectory  reentry 
point  into  the  layer  of 
branches. 

mm 

X'ri 

The  transformed  X7  coordi¬ 
nate  of  the  i*h  munition 
trajectory  reentry  point 
into  tho  layer  of  branches. 

mm 

Xsj 

The  X  ■  oordinate  of  an  end 
of  the  i**1  branch  after  rota¬ 
tion  tli rough  the  azimuth 
angle. 

mm 

Xt 

XTRAN 

The  distance  in  the  X  direc¬ 
tion  used  to  translate  the 
rotated  branch  coordinates 
into  the  first  octant  (see 

Table  1 J. 

mm 

Y'di 

YDET 

The  transformed  Y7  coordi¬ 
nate  of  the  detonation  for 
the  ilh  munition  trajectory. 

mm 
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SYMBOL  OR 
ABBREVIATION 

EQUIVALENT 

IN 

SIMULATION 

MODEL 

DEFINITION 

UNITS 

YEi 

YTRAJ(IOO) 

The  Y  coordinate  of  the  i^1 
munition  trajectory  entry 
point  into  the  layer  of 
branches. 

mm 

V'ci 

YTRAJP 

The  transformed  Y'  coordi¬ 
nate  of  the  i*h  munition 
trajectory  entry  point  into 
the  layer  of  branches. 

mm 

YJ 

BRANCH(6) 

The  Y  coordinate  of  an  end 
of  the  jth  branch. 

mm 

Y 1(1500) 
Y2(1500) 

The  transformed  Y' coordi¬ 
nate  of  an  end  of  the  jth 
branch. 

mm 

Y 

1  max 

YLIMIT 

The  Y  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment. 

mm 

YOi 

The  Y  coordinate  of  the  i*h 
munition  trajectory  exit, 
point  from  the  layer  of' 
branches. 

mm 

Yri 

The  Y  coordinate  of  the  i^*1 
munition  trajectory  reentry 
point  into  the  layer  of 
branches. 

mm 

Y*  ■ 

1  ri 

The  transformed  Y*  coordi¬ 
nate  of  the  ith  munition 
trajectory  reentry  point 
into  the  layer  of  branches. 

mm 

xiii 
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EQUIVALENT 

IN 

SYMBOL  OR  SIMULATION 

ABBREVIATION  MODEL _ DEFINITION _  UNITS 

Ygj  The  Y  coordinate  of  an  end  mm 

of  the  branch  after  rota¬ 
tion  through  the  azimuth 
angle. 

Yt  YTRAN  The  distance  in  the  Y  direc-  mm 

tion  used  to  translate  the 
r  otated  branch  coordinates 
into  the  first  octant  (see 
Table  f). 

YTRAJP  The  tr  ansformed  Y'  coordi-  mm 

nate  of  the  i*h  munition 
trajectory. 

Yjj  Yl(1500)  Tire  transformed  Y'  coordi-  mm 

nate  of  the  source  end  of  the 
jfh  branch. 

Yaj  Y2(1500)  The  transformed  Y1  coordi-  mm 

nate  of  the  terminal  end  of 
the  jth  branch. 

BOTTOM  The  Z  coordinate  of  the  bot-  mm 

tom  of  the  layer  of  branches. 

Z^  ZDET  The  transformed  z!  coordi-  mm 

nate  of  the  detonation  for 
the  ith  munition  trajectory. 

Zg^  ZTRAJ(lOO)  The  Z  coordinate  of  the  i^  mm 

munition  trajectory  entry 
point  into  the  layer1  of 
branches. 


xiv 
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(CONTINUED) 

EQUIVALENT 

IN 

SYMBOL  OR  SIMULATION 

ABBREVIATION  MODEL  DEFINITION  UNITS 


zEi 

ZTRAJP 

The  transformed  Z'  coordi¬ 
nate  of  the  i4*1  munition 
trajectory  entry  point  into 
the  layer  of  branches. 

mm 

Zh: 

ZHGT 

The  burst  height  for  the  i**1 
munition  trajectory. 

feet 

BRANCH(6) 

The  Z  coordinate  of  an  end 
of  the  jth  branch. 

mm 

z'j 

Zl(1500) 
Z2(l  500) 

The  transformed  z'  coordi¬ 
nate  of  an  end  of  the 
branch. 

mm 

f 

'max 

Z  LIMIT 

The  Z  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment. 

mm 

ZOi 

ZOUT 

The  Z  coordinate  of  the  i™ 
munition  trajectory  exit 
point  from  the  layer  of 
branches. 

mm 

zri 

The  Z  coordinate  of  the  i**1 
munition  trajectory  reentry 
point  into  the  layer  of 
branches. 

mm 

Z  ri 

The  transformed  z!  coordi¬ 
nate  of  the  i*h  munition 
trajectory  reentry  point 
into  the  layer  of  branches. 

mm 

xv 
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(CONCLUDED) 

EQUIVALENT 

IN 

SYMBOL  OR  SIMULATION 


ABBREVIATION 

MODEL 

DEFINITION 

UNITS 

Z*J 

The  Z  coordinate  of  an  end 
of  the  jth  branch  after  rota¬ 
tion  through  the  azimuth 
angle. 

mm 

z'ti 

ZTRAJP 

The  transformed  7J  coordi¬ 
nate  of  the  i*h  munition 
t  rajectory . 

mm 

zM 

ZK1500) 

The  transformed  z!  coordi¬ 
nate  of  the  source  end  of  the 
branch. 

mm 

Z'2J 

Z2(1500) 

The  transformed  Z'  coordi¬ 
nate  of  the  terminal  end  of 
the  jth  branch. 

mm 

0 

THETA 

The  azimuth  angle  for  the 
munition  trajectory  paths. 

degrees 

cp 

PHI 

The  elevation  angle  for  the 
munition  trajectory  paths. 

degrees 

\  \  i 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS  -  SIMULATION  MODE  I 


SYMBOL  OR 
ABBREVIATION 

EQUIVALENT 
IN  MATH 
MODEL 

DEFINITION 

UNITS 

A 

A 

The  distance  between  the 
branch  ends  along  the  trans¬ 
formed  Z7  axis. 

mm 

AVDIAMU500) 

An  array  containing  the 
average  diameters  of  the 
branches. 

mm 

AVRAD 

Rbj 

The  average  radius  of  the 
branch. 

mm 

AZIMUTH 

9 

The  azimuth  angle  for  the 
munition  trajectory  paths. 

radian? 

B 

B 

The  distance  between  the 
branch  ends  along  the  trans¬ 
formed  Y7  axis. 

mm 

BBDIAM 

Dd 

The  minimum  branch  dia¬ 
meter  required  for  detona¬ 
tion  by  the  munition  body. 

mm 

BDIAM 

2*m‘ 

The  diameter  of  the  muni¬ 
tion  body. 

mm 

BFDIAM 

The  minimum  branch  dia¬ 
meter  required  for  detona¬ 
tion  by  the  munition  fuze. 

mm 

BMT 

The  difference  in  height  (or 
the  difference  in  Z  coordi¬ 
nates)  between  the  top  and 
the  bottom  of  the  layer  of 
branches. 

mm 

BOTTOM 

zb 

The  Z  coordinate  of  the  bot¬ 
tom  of  the  layer  of  branches. 

mm 
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SYMBOL  OR 
ABBREVIATION 

EQUIVALENT 
IN  MATH 
MODEL 

DEFINITION 

UNITS 

BRAD 

Rm 

The  radius  of  the  munition 
body. 

mm 

BRANCH(6) 

XJ 

YJ 

Zj 

An  array  containing  the  X, 

Y,  and  Z  coordinates  of  the 
branch  ends. 

mm 

BRDIAM 

The  average  diameter  of  the 
branch  which  caused  the 
detonation. 

mm 

C 

C 

A  value  used  in  computing 
the  distance  in  the  trans¬ 
formed  Y'-Z'  plane  between 
the  munition  trajectory  and 
the  branch  centerline  for 

Cases  2  and  3. 

mm 

CFREQU50) 

Ck 

The  cumulative  frequency 
for  each  interval  of  the 
burst  height  distribution 
table. 

percent 

COSP 

CO  Sep 

The  cosine  of  the  elevation 
angle. 

none 

COST 

COS0 

The  cosine  of  the  azimuth 
angle. 

none 

CTCP 

cosGcoscp 

The  cosine  of  the  azimuth 
angle  multiplied  by  the  co¬ 
sine  of  the  elevation  angle. 

none 
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(CONTINUED) 


EQUIVALENT 
S\  M  HOI  .OR  IN  MATH 

ABBREVIATION  MODEL 

CTSP  cosSsincp 

D  D8 

DET(IOO) 

DIAM  2Rbj 

DUMMY 

DI 

FACTRB  (Rbj  +  Rm)s 

FACTRF  (Rbj  +  Rf)2 

FDIAM  2Rf 


DEFINITION 


The  cosine  of  the  azimuth 
angle  multiplied  by  the  sine 
of  the  elevation  angle. 

The  square  of  the  distance  in 
the  transformed  Y  - Z '  plane 
between  the  munition  trajec¬ 
tory  and  the  branch  center- 
line. 

An  array  containing  the  burst 
heights  for  100  munition 
trajectories. 

The  average  diameter  of  the 
branch. 

The  random  number  genera¬ 
tor  starting  value. 

An  intermediate  value  used  :.n 
establishing  the  burst  height 
intervals. 

The  square  of  the  sum  of  the 
average  radius  of  the  branch 
and  the  radius  of  the  muni¬ 
tion  body. 

The  square  of  the  sum  of  the 
average  radius  of  the  branch 
and  the  radius  of  the  muni¬ 
tion  fuze. 

The  diameter  of  the  muni¬ 
tion  fuze. 


UNITS 

none 


v__  2 

mm 


feet 


mm 


none 


none 


mm2 


mm2 


mm 
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SYMBOL  OR 
ABBREVIATION 

EQUIVALENT 
IN  MATH 
MODEL 

DEFINITION 

UNITS 

FRAD 

Rf 

The  radius  of  the  munition 
fuze. 

mm 

FREQ(150) 

Fk 

An  array  containing  the  num¬ 
ber  of  munition  detonations 
tor  each  burst  height  interval. 

none 

I 

An  internal  program  variable. 

none 

INC 

k 

The  sequential  burst  height 
interval  (i.  e.  ,  1  identifies 
ground  burst?,  2  identifies 
the  first  interval  above  the 
ground,  etc.  ). 

none 

I  PRINT 

The  print  option: 

0  =  print  the  burst  height 
and  branch  diameter  for 
each  munition  detona¬ 
tion  and  the  burst  height 
distribution  table. 

1  -  print  the  burst  height 
distribution  table. 

none 

IPUNCH 

• 

The  punch  option: 

0  =  do  not  punch  cards. 

1  =  punch  cards  for  the 
burst  height  intervals, 
the  number  of  detona¬ 
tions  in  each  interval, 
and  the  cumulative  fre¬ 
quency  for  each  interval. 

none 

ITHETA 

0/90  +  1 

An  integer  number  which  is  a 
function  of  the  azimuth  angle. 

none 

\  x 
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SYMBOL  OR  IN  MATH 
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IX 


J 

K 

KEEP 


KEOF 


KEPDET(300) 


KEPDIAM(300) 


KTGRD 

KTR 

L 


The  number  of  munition  deto-  none 

nations  for  the  burst  height 
interval. 

An  internal  program  variable.  none 

An  internal  program  variable.  none 

A  counter  for  the  number  of  none 


detonations  occurring  along 
each  munition  trajectory  in 
the  layer  of  branches. 

The  argument  to  System  Sub-  none 
routine  EOF  which  exercises 
control  over  the  program 
when  an  end-of-file  is  read. 

An  array  containing  the  burst  feet 
heights  (or  Z  coordinates) 
for  detonations  occurring 
along  egch  munition  trajec¬ 
tory  in  the  layer  of  branches. 

An  array  containing  the  branch  mm 
diameters  for  detonations 
occurring  along  each  munition 
trajectory  in  the  layer  of 
branches. 

The  number  of  ground  bursts.  none 
The  number  of  air  bursts.  none 

An  internal  program  variable.  none 


xxi 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS  -  SIMULATION  MODEL 

(CON TINIJ  ED) 


SYMBOL  OR 
ABBREVIATION 

EQUIVALENT 
IN  MATH 
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DEFINITION 

UNITS 

LENGTH 

(Aa  +  B3)^ 

The  length  of  the  branch  in 
the  transformed  x' -z! 
plane. 

mm 

LENSQR 

A3  *  B2 

The  square  of  the  length  of 
the  branch  in  the  transformed 
X'-Z.'  plane. 

mm2 

LEVEL 

A  seque  ntial  number  which 
identifies  the  layers  of 
branch  data  (i.  e.  ,  1  iden¬ 
tifies  da.ta  for  the  top  layer, 

2  identifies  data  for  the  next 
layer  from  the  top,  etc.  ). 

none 

LL 

An  internal  program  variable. 

none 

MIN 

A  subscript  used  to  identify 
the  branch  end  which  is 
closest  to  the  munition  tra¬ 
jectory  [see  R(MIN)]. 

none 

NBRAN 

The  number  of  branches  de¬ 
scribed  in  the  layer  on  the 
branch  input  data  tape. 

none 

NMUN 

The  munition  identification 
number. 

none 

NRNX 

An  intermediate  value  used  to 
determine  uniform  random 
integers  between  1  and  10. 

none 

NRNXY 

An  intermediate  value  used  to 
determine  uniform  random 
integers  between  1  and  10. 

none 
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SYMBOL  OR 
ABBREVIATION 

IN  MATH 
MODEL 

DEFINITION 

UNITS 

PHI 

CP 

The  elevation  angle  for  the 
munition  trajectory  paths. 

degrees 

PH  ISA  V 

A  save  location  for  the  ele¬ 
vation  angle. 

degrees 

PHISAVE 

A  save  location  for  the  ele¬ 
vation  angle. 

degrees 

PX 

The  X'  coordinate  of  the  deto¬ 
nation  measured  from  the 
ground  to  the  nose  of  the 
munition. 

mm 

R(2 ) 

Rai 

R3a 

An  array  containing  the 
squares  of  the  distances  in 
the  transformed  Y'-Z'  plane 
between  the  munition  trajec¬ 
tory  and  each  end  of  the 
branch. 

mm2 

REPEAT 

A  flag  which  when  set  to  1 
indicates  that  the  munition 
trajectory  must  be  reentered 
into  the  layer  of  branches. 

none 

R(MIN) 

p2 

^  min 

The  square  of  the  smallest  of 
the  distances  in  the  trans¬ 
formed  Y/-Z/  plane  between 
the  munition  trajectory  and 
each  end  of  the  branch. 

mm3 

RNX 

RXi 

A  random  integer  between  1 
and  10  selected  from  a  uni¬ 
form  distribution. 

none 
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DEFINITION 

RNY 

Ryi 

A  random  integer  between  1 
and  10  selected  from  a  uni¬ 
form  distribution. 

SINP 

sincp 

The  sine  of  the  el  'vation 
angle. 

SINT 

sin9 

The  sine  of  the  azimuth 
angle. 

SITE 

The  identifying  notation  for 
the  forest  environment. 

SPOCP 

tancp 

The  tangent  of  the  elevation 
angle. 

SSQ 

R2max  "  £)2 

The  difference  between  the 

square  of  the  largest  of  the 
distances  in  the  transformed 


Y'-Z'  plane  between  the  muni¬ 
tion  trajectory  and  each  end 
of  the  branch  and  the  square 
of  the  distance  in  the  trans¬ 
formed  Y'-Z1  plane  between 
the  munitiqn  trajectory  and 
the  branch  centerline. 

SS2  S22  The  square  of  the  distance  in 

the  transformed  Y'-Z'  plane 
along  the  centerline  of  the 
branch  and  between  one  of  the 
ends  of  the  branch  and  the 
point  at  which  the  normal  from 
the  munition  trajectory  inter¬ 
sects  the  branch  centerline. 


UNITS 

none 


none 


none 


none 


none 


mm2 


mm2 
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SYMBOL  OR 
ABBREVIATION 

STAT(8) 

STCP 

STSP 

SI 


S2 


THEDET 

THEDIAM 


EQUIVALENT 
IN  MATH 

MODEL  DEFINITION 


An  array  containing  statisti¬ 
cal  values  returned  by  Sys¬ 
tem  Subroutine  BDS. 

sinScoscp  The  sine  of  the  azimuth  angle 
multiplied  by  the  cosine  of 
the  elevation  angle. 

sinGsincp  The  sine  of  the  azimuth  angle 
multiplied  by  the  sine  of  the 
elevation  angle. 

Sx  The  distance  in  the  trans¬ 

formed  Y/-Z/  plane  along  the 
centerline  of  the  branch  and 
between  one  of  the  ends  of 
the  branch  and  the  point  at 
which  the  normal  from  the 
munition  trajectory  inter¬ 
sects  the  branch  centerline. 


S2  The  distance  in  the  trans¬ 

formed  Y/-Z/  plane  along  the 
centerline  of  the  branch  and 
between  one  of  the  ends  of 
the  branch  and  the  point  at 
which  the  normal  from  the 
munition  trajectory  inter¬ 
sects  the  branch  centerline. 


An  intermediate  location  for 
storing  the  burst  height  of 
the  munition  trajectory. 

The  branch  diameter  which 
caused  the  munition  detonation. 


UNITS 

various 

none 

none 

mm 


mm 


feet 

mm 
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SYMBOL  OR 
ABBREVIATION 

THETA 

TOP 

TRASH 

XCELL 

XDATA 

XDET 

XINT 

XLIMIT 

XLIMO 

XLIMOO 


EQUIVALENT 
IN  MATH 

MODEL  DEFINITION 


0  The  azimuth  angle  for  the 

munition  trajectory  paths. 

The  X  coordinate  of  the  top 
of  the  layer  of  branches. 

A  dummy  argument  for  the 
system  uniform  random  num¬ 
ber  generator. 

The  minimun  X  coordinate  of 
a  rectangle  which  subdivides 
the  top  of  the  forest  environ¬ 
ment. 

The  upper  boundary  of  the 
burst  height  interval. 

X^  The  transformed  X/  coordi¬ 

nate  of  the  detonation  for  the 
munition  trajectory. 

The  burst  height  interval. 

Xmax  The  X  coordinate  which  de¬ 

fines  a  boundary  of  the  forest 
environment. 

The  X  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment  divided  by  10 
(i.e.,  XLIMIT/ 10). 

The  X  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment  divided  by  100 
(i.e.  ,  XLIMIT /  100). 


UNITS 

degrees 

mm 

none 

mm 

feet 

mm 

feet 

mm 

mm 

mm 


x.wi 
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S\  MMOL  OK 
ABBREVIATION 

XNUM 

XOUT 

XPLUS 


XTRAJ(IOO) 

XT  RAN 

XTRANC 


EQUIVALENT 
IN  MATH 

MODEL  DEFINITION 


Tin*  total  number  of  munition 
trajectories  evaluated  for 
each  elevation  angle  (always 
400). 

XOi  The  X  coordinate  of  the  muni¬ 

tion  trajectory  exit  point  from 
the  layer  of  branches. 

The  distance  traveled  by  the 
munition  along  the  X  axis 
from  the  entry  point  into  the 
top  of  the  layer  to  the  exit 
point  through  the  plane  de¬ 
fining  the  bottom  of  the 
layer. 

Xgi  An  array  containing  the  X 

coordinates  of  the  munition 
trajectory  entry  points  into 
the  layer. 

Xt  The  distance  in  the  X  direc¬ 

tion  used  to  translate  the 
rotated  branch  coordinates 
into  the  first  octant  (see 
Table  I). 


X^coscp  The  distance  in  the  X  direc¬ 
tion  used  to  translate  the 
rotated  branch  coordinates 
into  the  first  octant  multi¬ 
plied  by  the  cosine  of  the 
elevation  angle. 


xxv  ii 


UNITS 

none 


mm 


mm 


mm 


mm 


mm 
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UNITS 

XTRANS 

Xtsincp 

The  distance  in  the  X  direc¬ 
tion  used  to  translate  the 
rotated  branch  coordinates 
into  the  first  octant  multi¬ 
plied  by  the  sine  of  the 
elevation  angle. 

mm 

XX 

The  lower  boundary  of  the 
burst  height  interval. 

feet 

Xl(1500) 

An  array  containing  the  trans¬ 
formed  X'  coordinates  for  the 
source  ends  of  the  branches. 

mm 

X2(1500) 

X3j 

An  array  containing  the  trans¬ 
formed  X1  coordinates  for  the 
terminal  ends  of  the  branches. 

mm 

YCELL 

The  minimum  Y  coordinate  of 
a  rectangle  which  subdivides 
the  top  of  the  forest  environ¬ 
ment. 

miu 

YDET 

Y'di 

i 

The  transformed  Y*  coordi¬ 
nate  of  the  detonation  for  the 
munition  trajectory. 

mm 

YLIMIT 

Y 

1  max 

The  Y  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment. 

mm 

YLIMO 

The  Y  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment  divided  by  10 
(i.  e. ,  YLIMIT / 10). 

mm 
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UNITS 

YLIMOO 

The  Y  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment  dividdd  by  100 
(i.e.,  YLIMIT/100). 

mm 

YTRAJ(IOO) 

YEi 

An  array  containing  the  Y 
coordinates  of  the  munition 
trajectory  entry  points  into 
the  layer  of  branches. 

mm 

YTRAJP 

YEi 

The  transformed  Y/  coordi¬ 
nate  of  the  munition  trajec¬ 
tory  entry  point  into  the 
layer  of  branches. 

mm 

YTRAN 

Yt 

The  distance  in  the  Y  direc¬ 
tion  used  to  translate  the 
rotated  branch  coordinates 
into  the  first  octant  (see 

Table  I). 

mm 

YH1500) 

v'.i 

An  array  containing  the  trans¬ 
formed  Y;  coordinates  of  the 
source  ends  of  the  branches. 

mm 

Y2U500) 

Y2J 

An  array  containing  the  trans¬ 
formed  Y/  coordinates  of  the 
terminal  ends  of  the  branches. 

mm 

ZDET 

z'di 

The  transformed  z!  coordi¬ 
nate  of  the  detonation  for  the 
munition  trajectory. 

mm 

ZHGT 

zhi 

The  burst  height  for  the  muni¬ 
tion  trajectory, 

feet 
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SYMBOL  OR 
ABBREVIATION 

ZLIMIT 

ZOUT 

ZTRAJ(IOO) 

ZTRAJP 

Zl(1500) 

Z2(1500) 


EQUIVALENT 
IN  MATH 

MODEL  DEFINITION 


7 

max 


Z 


Oi 


zEi 


The  Z  coordinate  which  de¬ 
fines  a  boundary  of  the  forest 
environment. 

The  Z  coordinate  of  the  muni¬ 
tion  trajectory  exit  point  from 
the  layer  of  branches. 

An  array  containing  the  Z 
coordinates  of  the  munition 
trajectory  entry  points  into 
the  layer  of  branches. 

The  transformed  z!  coordi¬ 
nate  of  the  munition  trajec¬ 
tory  entry  point  into  the 
layer  of  branches. 

An  array  containing  the  trans¬ 
formed  7.!  coordinates  of  the 
source  ends  of  the  branches. 

An  array  containing  the  tran- 
formed  Z‘  coordinates  of  the 
terminal  ends  of  the  branches. 


w\ 


UNITS 

mm 


mm 


mm 


mm 


mm 


mm 


SECTION  I 


INTRODUCTION 


The  Burst  Height  Distribution  (BHD)  Program  was  designed  to 
compute  and  display  burst  height  distributions  for  munitions  aerially 
delivered  into  forest  environments.  The  computer  program  was  written 
in  FORTRAN  IV  and  is  presently  operational  on  the  Control  Data 
Corporation  6600  computer  system  at  Eglin  Air  Force  Base,  Florida. 

The  program  uses  as  input  the  source  and  terminal  X,  Y,  and  Z 
coordinates  and  the  average  diameters  of  branches  surveyed  at  actual 
forested  sites.  Because  of  the  large  number  of  branches  included  in  the 
characterization  of  the  forest  environments,  all  of  the  branch  coordi¬ 
nates  and  diameters  cannot  be  read  into  the  computer  central  processor 
at  one  time.  An  input  tape  must  be  prepared  with  the  branch  data  seg¬ 
mented  into  layers.  The  branch  data  for  the  top  layer  are  then  read  into 
the  computer  central  processor,  and  the  munition  trajectories  are 
examined  for  possible  encounters  with  the  branches.  For  those  munition 
trajectories  which  pass  through  the  top  layer  and  do  not  result  in  a 
munition  detonation,  the  coordinates  of  the  exit  point  from  the  bottom  or 
side  of  the  layer  are  determined.  If  the  munition  trajectory  exits  the 
side  of  the  layer,  the  trajectory  is  reentered  into  the  layer  at  the  same 
height.  If  the  trajectory  exits  the  bottom  of  the  layer,  the  exit  coordi¬ 
nates  from  the  top  layer  are  used  as  entry  coordinates  for  the  second 
layer.  This  process  is  repeated  until  all  layers  have  been  examined  or 
until  the  munition  detonates.  The  procedure  permits  determination  of 
the  munition/ branch  encounters  for  all  of  the  munition  trajectories  by 
inputting  the  branch  data  into  the  central  processor  only  one  time. 

The  munitions  travel  along  straight-line  trajectories  which  are 
randomly  selected.  If  an  encounter  between  a  branch  and  the  munition 
body  occurs,  the  program  determines  if  the  branch  is  large  enough  to 
cause  the  munition  to  detonate.  If  the  branch  is  large  enough,  the  burst 
height  is  calculated.  If  the  branch  is  too  small,  the  program  next  deter¬ 
mines  if  the  munition  fuze  encounters  the  branch.  If  the  fuze  encounters 
the  branch,  the  burst  height  is  calculated,  and  if  the  fuze  does  not  en¬ 
counter  the  branch,  the  trajectory  is  continued.  After  100  trajectories 
are  examined  from  each  of  four  azimuth  angles  (i.  e.,  0,  90,  180,  and 
270  degrees)  and  from  the  elevation  angle  specified  as  input,  the  cumu¬ 
lative  burst  height  distribution  for  the  munition  is  computed,  printed, 
and  optionally  punched  as  output. 
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This  volume  contains: 


•  A  detailed  description  of  the  mathematical  model  which  was 
used  as  a  basis  in  developing  the  program. 

•  A  complete  set  of  flowcharts  depicting  the  logic  used  in  the 
program. 

•  A  complete  FORTRAN  IV  source  listing  of  the  program. 

•  A  detailed  description  of  the  simulation  coding  which  was 
employed  in  developing  the  program. 

•  Listings  of  the  symbols  and  abbreviations  which  are  used  in 
the  mathematical  model  and  the  computer  program. 

Detailed  descriptions  of  the  input  variables  required  to  properly 
execute  the  program,  instructions  for  placing  the  input  variables  on 
punch  cards  and  for  arranging  the  punch  cards  in  proper  order,  descrip¬ 
tions  and  definitions  of  the  output  available  from  the  program,  and  a 
description  of  a  sample  case  are  contained  in  Volume  I,  the  User  Manual, 
of  this  report. 
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SECTION  II 


MATHEMATICAL  MODEL 


The  mathematical  model  is  divided  into  five  principal  parts.  The 
first  of  these  determines  the  entry  coordinates  for  each  of  the  m  .nition 
trajectories  that  are  to  be  examined  as  they  pass  through  the  top  of  the 
forest  environment.  The  second  part  translates  and  rotates  the  coordi¬ 
nates  of  the  branches  so  that  they  are  in  the  first  octant  and  so  that  the 
new  x'  axis  is  parallel  to  the  path  of  the  munition  trajectory.  The  third 
part  determines  subsequent  munition  trajectory  exit  and  reentry  points 
as  they  pass  through  the  layers  of  the  forest  environment.  The  fourth 
part  calculates  the  distances  between  each  of  the  munition  trajectories 
and  branch  centerlines  and  determines  the  munition/branch  encounters 
with  the  maximum  height  that  will  detonate  a  munition  traveling  along 
each  of  the  trajectory  paths.  The  fifth  part  prepares  cumulative  burst 
height  distribution  tables  for  the  400  munition  trajectories  delivered  from 
the  four  azimuth  angles  (i.e.,  100  each  from  0,  90,  180,  and  270  degrees) 
and  from  a  specified  elevation  angle. 


DETERMINATION  OF  INITIAL  TRAJECTORY  ENTRY  COORDINATES 
INTO  TOP  LAYER 

The  source  and  terminal  X,  Y,  and  Z  coordinates  of  the  branches 
in  the  forest  environment  are  originally  specified  in  the  first  octant  of 
a  cartesian  coordinate  system,  and  the  limits  of  the  environments  are 
defined  by  Xmax,  Ymax,  and  Zmax,  as  shown  in  Figure  1.  The  branch 
coordinates  and  average  diameters  are  contained  6n  magnetic  tape,  and 
the  tape  is  subdivided  into  layers  with  data  for  up  to  1,  500  branches  con¬ 
tained  in  each  layer.  The  entry  coordinates  for  each  of  the  munition 
trajectories  which  are  to  be  examined  from  a  given  azimuth  angle  and 
elevation  angle  are  determined  by  subdividing  the  X-Y  plane  of  the  forest 
environment  into  100  rectangles  and  randomly  selecting  entry  coordi¬ 
nates  which  will  pass  through  each  of  the  rectangles.  The  expressions 
used  in  the  model  are: 


(p-D(Xmax)  RXi(Xmax) 
XEi  =  10  +  100 


(1) 
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Figure  1.  Pictorial  Representation  of  Layers  Used  in  Forest 
Environment  and  Trajectory  Paths 
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,,  (q-l)(Ymax>  ^Yi^max^ 
YEi  =  - - + 


10 


100 


Ze  i  s  Zmax 


(2) 

(3) 


w>, ere  p  is  tlio  sequential  number  of  the  rectangle  along  the  X  axis,  q  is 
the  sequential  number  of  the  rectangle  along  the  Y  axis,  Xgj,  Ygi.  and 
ZEi  are  the  X,  Y,  and  Z  coordinates  of  the  i*h  munition  trajectory  entry 
point,  Xmax,  Ymax,  and  Zmax  are  the  X,  Y,  and  Z  boundaries  of  the 
forest  environment  (see  Figure  1),  and  and  Ry[  are  random  integers 
between  1  and  10  selected  from  a  uniform  distribution. 


TRANSLATION  AND  ROTATION  OF  BRANCH  COORDINATES 

A  time  saving  feature  of  the  BHD  Program  involves  a  transforma¬ 
tion  of  the  original  coordinate  axes  to  permit  determination  of  the 
munition/branch  encounters  in  two  rather  than  three  dimensions.  The 
transformation  is  made  so  that  the  new  X7  axis  is  parallel  to  the  path  of 
the  munition  trajectory,  and  all  subsequent  munition/branch  encounter 
calculations  are  performed  using  branch  coordinates  projected  onto  the 
new  Y'-Z'  plane.  The  equations  for  the  first  rotation  about  the  Z  axis 
through  the  azimuth  angle  9  and  for  translation  of  the  coordinates  into 
the  first  octant  are  (see  Figure  2): 


Xsj  =  XjCosB  +  YjSinS  +  X* 

Ygj  =  -XjsinS  +  YjCosB  +  Yt 
Zsj  =  Zj 

where  9  is  as  defined  above,  Xj,  Yj,  and  Zj  are  the  coordinates  of  the 
jth  point  relative  to  the  original  X,  Y,  and  Z  axes,  X0j,  YSj,  and  ZSj  are 
the  new  translated  and  rotated  coordinates,  and  X^  and  Y^  are  the  X  and 
Y  values  used  to  translate  the  new  coordinates  into  the  first  octant.  In 
matrix  form,  this  first  transformation  is  given  by! 
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TRANSFORMED 
POINT  (XsjYsj) 


ORIGINAL 
POINT  (XjYj) 


YjCosQ 


F 


xsj" 

COS0 

sinO 

0 

V 

xt 

Ysj 

= 

-sinfl 

cost) 

0 

Y.i 

+ 

Yt 

ZsJ_ 

0 

0 

zb 

0 

The  second  rotation  is  about  the  newly  established  Ys  axis  through 
the  elevation  angle  cp,  which  establishes  the  X',  Y',  and  Z*  coordinate 
axes.  The  equations  for  this  second  rotation  are  (see  Figure  3): 


XXj  =  XSjCoscp  -t  Zgjsincp 


Y  =  Y 
J  s.i 


Z7j  =  -XS|Sincp  +•  Zs)coscp 


In  matrix  form,  this  second  transformation  is  given  by: 


fx/n 

XJ 

r 

coscp 

0 

sincp 

V 

Y'j 

r 

0 

1 

0 

YSJ 

Z/J 
L  JJ 

-sincp 

0 

coscp 

Z  o  i 

_ 

The  composite  equations  encompassing  both  rotations  and  the  transla¬ 
tion  are  expressed  in  matrix  form  as: 


coscp 

0 

sincp 

COS  0 

sinP 

(f 

X 

= 

0 

1 

0 

-sin0 

COS0 

0 

YJ 

.ZL 

-sincp 

0 

coscp 

J 

0 

tm 

0 

1 

ZJ 

- 

coscp  0  sincp 

X 

+ 

0  1  0 

Yt 

-sincp  0  coscp 

0 
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Multiplying  the  two  square  matrices  results  in  the  following: 


a 

cos0coscp  sin8cos(p  sincp 

A 

coscp  0  sincp 

Xt‘ 

u 

r 

-sin0  cos0  0 

yj 

+ 

0  i  0 

Yt 

A 

-cosOsincp  -sinSsincp  coscp 

_ 

A 

-sincp  0  coscp 
— 

0 

_  _ 

Thus,  the  equations  expressing  the  newly  established  x',  Y',  and 
Z'  coordinates  as  functions  of  the  original  X,  Y,  and  Z  coordinates,  the 
azimuth  angle  0,  and  the  elevation  angle  cp  are: 


x'j  =  XjCos8coscp  +  Yjsin6coscp  +  Zjsincp  +  Xtcoscp  (4) 

Y#j  =  -XjSinfl  +  YjCosQ  +  Yt  (5) 

=  *XjCos8sincp  -  YjSinOsincp  +  Zjcoscp  -  Xj-sincp  (6) 

where  all  variables  are  as  defined  above.  The  values  of  Xt  and  Yt  used 
in  the  program  vary  with  the  azimuth  angle  and  are  presented  in  Table  I. 


TABLE  I.  TRANSLATION  VALUES  Xt  AND  Yt 
USED  IN  BHD  PROGRAM 


AZIMUTH  ANGLE 
(degrees) 

Xt 

Yt 

0 

0 

0 

90 

0 

Y  max 

180 

V 

Amax 

Ymax 

270 

^max 

0 
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DETERMINATION  OF  SUBSEQUENT  TRAJECTORY  EXIT  AND 
REENTRY  COORDINATES 


* 


The  X,  Y,  and  Z  coordinates  for  each  of  the  munition  trajectories 
as  they  exit  the  layers  of  branches  through  the  bottom  or  side  are  calcu¬ 
lated  using  the  expressions  (see  Figure  4): 


For  Munitions  Exiting  the 

For  Munitions  Exiting  the 

Layer  of  Branches  at  the 

Layer  of  Branches  at  the 

Bottom: 

Side: 

(ZEi  "  Zh) 

Xoi  =  +  tancp 

xOi  =  ^max 

(7) 

YOi  =  YEi 

YOi  =  YEi 

(8) 

N 

O 

II 

N 

cr 

zOi  =  zEi  ‘  (Xmax  '  xEi)tan(P 

(9) 

where  Xmax»  XEi,  YEp  and  ZEi  are  as  defined  above,  Xoi,  Yqp  and 
Zoi  are  the  coordinates  of  the  munition  trajectory  as  it  exits  the  layer, 

Zb  is  the  Z  coordinate  of  the  bottom  of  the  layer,  and  cp  is  the  elevation 
angle. 

If  the  munition  trajectory  exits  the  layer  through  the  side  rather 
than  through  the  bottom,  it  must  be  reentered  into  the  layer.  The  reentry 
process  is  repeated  until  the  trajectory  exits  the  bottom  of  the  layer,  and 
the  reentry  coordinates  are  determined  using  the  expressions. 


Xri  =  0.  0 

(10) 

1, 

Yri  =  ^N^Ymax^ 

(11) 

i? 

F 

t 

\ 

i 

zri  =  zOi 

(12) 

F 

I 

where  Ymax  and  Zqi  are  as  defined  above,  -Xri,  Yr[ 
reentry  coordinates  for  the  i*h  munition  trajectory, 
number  selected  from  a  uniform  distribution  with  0. 

,  and  Zr[  are  the 
and  Rn  is  a  random 

0  <  RN  <  1.  0. 
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Figure  4.  Relationship  of  Entry  and  Exit  Coordinates  as  Trajectorie 

Pass  Through  Layers  of  Branches 


When  the  munition  trajectory  exits  the  bottom  of  the  layer,  the 
reentry  coordinates  ore  determined  using  the  expressions: 


xri  =  xQi 

(13) 

Yri=  YEi 

(14) 

^ri  *  ^Oi 

(15) 

where  all  of  the  variables  are  as  defined  above. 

For  the  top  layer  of  branches,  the  entry  point  coordinates  of  each 
mm.  tion  trajectory  are  rotated  through  the  elevation  angle  cp  and  around 
t’.e  V  axis.  The  expressions  used  to  perform  this  rotation  are  (see  Fig- 

1  ' /: 


Xgi  =  Xgjcoscpf  ZgjSinqp 

(16) 

Y'ei  =  YEi 

(17) 

ZEi  =  -XEisincp  +  ZEiCoscp 

(18) 

where  XEi,  YEi,  and  ZEi  are  the  entry  coordinates  for  the  i^  munition 
trajectory  after  rotation,  and  all  other  variables  are  as  defined  above. 
For  each  succeeding  layer  of  branches  and  for  any  munition  trajectory 
which  must  reenter  a  layer,  the  reentry  coordinates  are  rotated  using 
the  expressions: 

XVi  =  Xr[Coscp  +  ZrjSincp 

(19) 

Y'  •  =  Y 

1  ri  1  ri 

(20) 

Z'vi  =  -X,.jsincp  +  ZrjCoscp 

(21) 

where  xVi.  Y^,  an(*  ^Vi  are  the  reentry  coordinates  for  the  i**1  munition 
trajectory  after  rotation,  and  all  other  variables  are  as  defined  above. 
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Figure  5.  Relationship  of  Original  and  Transformed  Munition 
Trajectory  Entry  Point  Coordinates  After  Rotation  in 
X-Y  Plane  Through  Elevation  Angle  cp 


DETERMINATION  OF  MUNITION /BRANCH  ENCOUNTERS  AND  BURST 
HEIGHTS 

The  munition/branch  encounters  are  evaluated  layer  by  layer  until 
each  munition  has  detonated  or  has  exited  the  bottom  layer  of  the  forest 
environment  without  detonating.  A  layer  of  branch  input  data  is  read 
into  core,  and  the  400  munition  trajectories  (100  from  each  of  the  four 
azimuth  angles)  are  evaluated  against  each  of  the  branches  in  the  layer. 
The  munition/branch  encounters  are  determined  by  considering  one  of 
the  following  three  cases: 


Case  1 


For  this  case,  the  branch  is  parallel  to  the  munition  trajectory 
(i.e.  ,  the  source  and  terminal  Y;  and  Z'  coordinates  for  the  branch  are 
identical,  as  shown  in  Figure  6).  The  square  of  the  distance  between  the 
munition  trajectory  and  the  branch  centerline  is  calculated  using  the 
expression: 


D2  =  (Y'j  -  Yxi)2  +  (Zj  -  z'Ti)3 


(22) 


where  I)  is  the  distance  between  the  munition  trajectory  and  the  branch 
centerline,  Y7j  and  Z7j  are  the  translated  coordinates  of  either  end  of  the 
jth  branch,  and  Y^i  and  Z  are  the  translated  entry  or  reentry  coordi¬ 
nates  of  the  i*h  munition  trajectory.  (In  this  and  subsequent  munition/ 
branch  encounter  computations,  comparisons  are  made  using  squared 
quantities.  This  technique  has  been  employed  to  reduce  the  computer 
running  time  for  the  program.  )  For  this  case,  a  munition  body/brjnch 
encounter  occurs  if: 


D3<(Rbj+Rm)3  (23) 

and  a  munition  fuze/branch  encounter  occurs  if: 

D8  <  (Rbj  +  Rf)2  (24) 
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Figure  6.  Munition/ Branch  Encounter  -  Case  1 
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where  O  is  as  defined  above,  R^j  is  the  average  radius  of  the  branch, 
Rm  is  the  munition  body  radius,  and  Rf  is  the  munition  fuze  radius. 

Only  those  branches  with  an  average  radius  R^j  greater  than  the 
minimum  branch  radius  required  for  detonation  by  a  munition  fuze  are 
read  m  from  the  branch  input  tape.  Therefore,  a  detonation  caused  by 
a  munition  fuze/ branch  encounter  will  occur  if  Inequality  24  is  satisfied. 
In  addition,  a  detonation  caused  by  a  munition  body /branch  encounter 
will  occur  if  Inequality  23  is  satisfied  and  if: 


2Rbj  >  Dd 


(25) 


where  Rbj  is  as  defined  above,  and  Dd  is  the  minimum  branch  diameter 
required  for  detonation  by  a  munition  body. 

If  Inequalities  23  and  25  or  Inequality  24  are  satisfied,  a  munition 
detonation  is  recorded,  and  the  transformed  coordinates  for  the  detona¬ 
tion  are  determined  using  the  expressions: 


X'di  =  X'min 

(26) 

Y'di  =  Y'ti 

(27) 

Z'di  =  Z'xi 

(28) 

where  Y'xi  and  Z'xi  are  as  defined  above,  X'di,  Y'di,  and  Z'di  are  the 
transformed  coordinates  of  the  detonation  for  the  i*h  munition  trajectory, 

!  and  X'min  is  the  smaller  of  the  transformed  X'  coordinates  of  the  branch 

ends. 

|  The  burst  height  for  the  detonation  is  then  computed  using  the 

expression: 

Zhi  s  0.  003281  (X/diSincp+  Z'djcoscp)  (29) 

t 

y 
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where  x'di,  Z'^,  and  cp  are  as  defined  above,  and  Z^i  is  the  burst  height 
of  the  i*h  munition  detonation. 


Case  2 


For  this  case,  the  branch  is  not  parallel  to  the  munition  trajectory, 
and  the  normal  from  the  munition  trajectory  to  the  centerline  of  the  branch 
intersects  the  branch  between  its  ends  or  nodes  (see  Figure  7).  To 
satisfy  the  conditions  for  this  case,  the  following  relationship  must  exist: 


(Ramax  -  D8)  <  (A3  +  Ba)  (30) 

where  R3max  is  the  greater  of  Rax  and  R33 : 

R2  =  (Y^j  -  YTi)3  +  (Z\j  -  Z'Ti)s  (31) 

Ra2  =  (Y'sj  -  Y'Ti)2  +  (Zaj  -  Z'Ti)3  (32) 

A  =  Z\j  -  Zaj  (33) 

B  =  Yaj  -  Y#xj  |  (34) 


Y  ij  and  Z  are  the  transformed  Y/  and  z'  coordinates  of  one  branch  end, 
Yaj  and  Z2j  are  the  transformed  Y/  and  z'  coordinates  of  the  other  branch 
end,  D  is  the  distance  between  the  munition  trajectory  and  the  branch 
centerline,  D2  can  be  determined  using  the  equation  for  computing  the 
distance  between  a  point  and  a  line,  as  follows: 


(A  Y^i  +  BZ'Tl  +  C)2 

D  A8  +  B2 

(35) 

where: 

C  =  -Zxj(B)  -  Y'xj(A) 

(36) 

and  all  other  variables  are  as  defined  above.  If  Inequalities  23  and  25  or 
Inequality  24  are  satisfied,  a  munition  detonation  is  recorded,  and  the 
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transformed  coordinates  for  the  detonation  are  determined  using  the 
expressions: 


where: 


X'di 

Si  x'sj  +  .S2XXJ 

Sx  +  S2 

(37) 

Ydi 

II 

(38) 

Z'di 

s 

KJ 

II 

(39) 

Si 

=  /A2  +  B2  -  S3 

(40) 

S3 

=  -  D2 

(41) 

and  all  other  variables  are  as  defined  above.  The  burst  height  for  the 
detonation  is  then  computed  using  Equation  29. 


Case  3 


For  this  case,  the  branch  is  not  parallel  to  the  munition  trajectory, 
and  the  normal  from  the  trajectory  to  the  centerline  of  the  branch  does  not 
intersect  the  branch  between  its  ends  or  nodes  (see  Figure  8).  To  satisfy 
the  conditions  for  this  case,  the  following  relationship  must  exist: 


(R2max  '  1)2 '>  -  (A2  +  Bs) 


(42) 


where  all  variables  are  as  defined  above.  A  munition  body/branch 
encounter  occurs  if: 


R  min  <  <Rbj  +  Rm>2  (43) 

» 

and  a  munition  fuze/branch  encounter  occurs  if: 

RSmin  <  (Rbj  +  Rf)  (44) 


19 


2 


r 


wlici  e  li  mm  is  the  smaller  of  If  ,  and  K  and  R[>|,  Km,  and  K|  arc  as 
defined  above.  If  Inequalities  25  anti  48  or  Inequality  44  are  satisfied,  a 
munition  detonation  is  recorded,  and  the  transformed  coordinates  for  the 
detonation  arc  determined  using  the  expressions: 


/ 

xVi 

if  R\  <  Ras 

(45) 

X  di  ■ 

x'aj 

if  R22  < 

Y'di  =  Y  Ti 

(46) 

zdi  =  zTi 

(47) 

where  all  variables  are  as  defined  above.  The  burst  height  for  the  deto¬ 
nation  is  then  computed  using  Equation  29. 


COMPUTATION  OF  CUMULATIVE  BURST  HEIGHT  DISTRIBUTIONS 

After  the  burst  heights  for  the  400  trajectories  have  been  determined, 
the?  model  groups  the  values  into  intervals  as  specified  in  the  input  and 
determines  the  number  of  detonations  and  the  cumulative  frequency  for 
each  burst  height  interval.  The  expression  which  determines  the  cumu¬ 
lative  frequency  is: 


lOOFk 

Ck  -  C  k-1  +  400 


(48) 


where  is  the  cumulative  frequency  (in  percent)  for  the  k^h  interval, 
and  F^  is  the  number  of  detonations  in  the  k*h  interval. 
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SECTION  III 


FLOWCHARTS 


This  section  contains  the  two  sets  of  flowcharts  which  depict  the 
logical  structure  of  the  main  program  and  its  subroutine.  The  flowcharts 
are  based  upon  the  logical  intent  of  the  model  rather  than  on  displaying 
the  methods  used  to  code  the  specific  program  routines.  Program  MAIN 
(Figure  9)  computes  the  munition  burst  heights  for  each  of  the  400  muni¬ 
tion  trajectories,  and  Subroutine  CFRQ  (Figure  10)  computes,  prints, 
and  optionally  punches  the  cumulative  burst  height  distribution  for  the 
munition  and  elevation  angle  being  considered. 
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WRITE  HEADING  AND 
CASE  IDENTIFYING 
INFORMATION 


Figure  9.  Flowchart,  Program  MAIN 
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Figure  9.  (Continued) 
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Figure  9.  (Continued) 


Figure  9.  (Continued) 
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Figure  !).  (Continued) 
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Figure  9.  (Continued) 
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Figure  9.  (Concluded) 
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Figure  10.  Flowchart,  Subroutine  CFRQ 
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Figure  10.  (Concluded) 


Auction  iv 


SOURCE  LISTING 


Figure  11  presents  a  complete  source  listing  of  the  BHD  Program, 
which  was  designed  for  operation  on  the  Control  Data  Corporation  6600 
computer  system  at  Eglin  Air  Force  Base,  Florida. 
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PROGRAM  P7821  ( INPUT,OUTPUT,  TAPE 5= I NPUT, TAPE6=0UTPUT , TAPE1 5,  TAPE  10) 

THIS  PROGRAM  DETERMINES  BURST  HEIGHT  DISTRIBUTIONS  FOR 
MUNITIONS  AERIALLY  DELIVERED  INTO  FOREST  ENVIRONMENTS 

DIMENSION  XTRAJU00),YTRAJO00),ZTRAJ(100),BRANCH(6),R(2),Xia500) 
1,  Yl(  1500),  Zl(  1500),  X2(  1500),  Y2(  1500),  Z2(  1500),  AVDIAMC 1500),  DET(  100 
2),KEPDET(300),KEPDIAM(  300) 

DIMENSION  STATC8) 

COMMON  NMUN,  S I TE,  PHI  SAV,  TFETA,  X INT,  I  PUNCH,  I  PRINT,  DET 
COMMON  XI, Yl,  Z1,X2, Y2,Z2, AVDIAM 
INTEGER  REPEAT 

REAL  LENSQR,  LENGTH, KEPDET,KEPDI AM 
CALL  EOF(KEOF) 

8  ASSIGN  300  TO  KEOF 

READ  IN  ELEVATION  ANGLE,  X,  Y,  AND  Z  ENVIRONMENT  BOUNDARIES, 

RANDOM  NUMBER  SEED,  BURST  HEIGHT  INTERVAL,  PUNCH  OPTION  (0=DON’T 
PUNCH  CARDS,  1=PUNCH  CARDS),  PRINT  OPTION  (0=PRINT  DATA  FOR  EACH 
DETONATION  AND  BURST  HEIGHT  DISTRIBUTION  TABLE,  1=PRINT  BURST 
HEIGHT  DISTRIBUTION  TABLE),  MUNITION  BODY  DIAMETER,  MUNITION 
FUZE  DIAMETER,  MINIMUM  BRANCH  DIAMETER  REQUIRED  FOR  DETONATION 
BY  MUNITION  BODY,  MINIMUM  BRANCH  DIAMETER  REQUIRED  FOR  DETONATION 
BY  MUNITION  FUZE,  MUNITION  NUMBER,  AND  SITE  DESIGNATION 

READC5,  1000)  PHI,XLIMIT,  YLIMIT,  ZLIMIT,  DUMMY,  XINT,  I  PUNCH,  I  PRINT, 
1BDIAM,  FDIAM,BBDIAM,BFDIAM,NMUN,  SITE 

PROVIDE  TRANSFERS,  INITIALIZE  VARIABLES  AND  ARRAYS,  AND  PERFORM 
INITIAL  CALCULATIONS. 

PH I S AVE=PH I 
THETA=0. 

6  REWIND  15 
PHI=PHISAVE 
DO  1  1=1,100 

1  DET(I )=0 . 

DO  2  1  =  1,300 
KEPDIAM(I)=0. 

2  KEPDET(I)=0. 

ASSIGN  110  TO  KEOF 
KEEP=0 

IFOHETA.NE.O.  .AND.  IPRINT.EQ.  1)  GO  TO  7 

WR I  TE  C  6 , 2  0  0  0  )  S I  TE ,  NMUN ,  BD I  AM ,  PH  I ,  FD I  AM,  THETA,  DUMMY  ,  BBD I  AM, 

1BFDIAM 


Figure  11.  BHD  Program  Source  Listing 
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I F( IPRINT.EQ.  1)  GO  TO  7 
WRITE  (6,2010) 

7  KTR=0 
PHI  SAV  =  PHI 
PHI =-PHI 
XTRAN=0. 

YTRAN=0. 

I  THETA=  THETA/90.  +  1. 

GO  TO(  5, 3, 3, 4),  I  THETA 

3  IFCITHETA.EQ.3)  XTRAN=XLfMIT 
YTRAN=YLIMIT 

GO  TO  5 

4  XTRAN=XLIMIT 

5  FRAD=FDIAM/2. 

Rn\D-BDIAM/2. 

rz  IMOTH= .  0 1 745  329  5 !{THETA 
PHI=.017453295!!PHI 
S I NT=S INC AZIMUTH) 

SINP=SIN(PHI) 

COST=COS(  AZIMUTH) 

COSP=COS(PHI) 

CTCP=COST:£COSP 
CTSP=COST«SINP 
STSP=SINT5£SINP 
ST  CP-S I NT”COSP 
XTRANC=XTRAN:£COSP 
XTRANS=XTRANJ£SINP 
SPOCP=S I NP/COSP 
C 

C  INITIALIZE  RANDOM  NUMBER  GENERATOR,  SUBDIVIDE  X-Y  PLANE  INTO 
C  100  RECTANGLES,  AND  RANDOMLY  SELECT  TRAJECTORIES  PASSING 
C  THROUGH  EACH  RECTANGLE. 

C 

XLIM0=XLIMIT/10. 

YLIMO=YLIMIT/10. 

XLIMOO=XLIMIT/100. 

YLIM0O=YLIMIT/100. 

CALL  RANSETC DUMMY) 

DO  10  1=1,10 
XCELL=XL  IM0:£(  I  - 1) 

K=(I-1)5£10 
DO  10  J=l,  10 
YCELL=YLIM0}£(J-1) 

NRNXY=RANF(TRASH)»100 . 

NRNX=NRNXY/ 10 
RNX=NRNX 


Figure  11.  (Continued) 
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RNY=NRNXY-(NRNXS!10) 

L=K+J 

XTRAJ(L)=XCELL+RNX"XLIMOO 
YT  RAJ  C  L  )  =Y  CE  LL+RN  Y5CYL I  MO  0 
10  ZTRAJ(L)=ZLIMIT 

READ  BRANCH  DATA  INTO  CORE,  ELIMINATE  BRANCHES  WITH  DIAMETER  LESS 
THAN  MINIMUM  REQUIRED  FOR  DETONATION  BY  MUNITION  FUZE,  AND 
TRANSFORM  COORDINATES  OF  BRANCH  ENDS  TO  X'-Y'-Z'  COORDINATE 
SYSTEM  WITH  X'  AXIS  PARALLEL  TO  TRAJECTORY  PATHS. 

TOP=ZLIMIT 

15  READC  15)  LEVEL,NBRAN,  BOTTOM, BRANCH,  DI AM 
K=0 

BMT =BOTT  OM-TOP 
DO  30  1=1, NB RAN 
IF(DIAM.LT.BFDIAM)GO  TO  20 
K=K  +  1 

Xl(K)=BRANCH(l):cCTCP+BRANCH(2)“STCP+BRANCH(3)”SINP+XTRANC 
X2(K)=B  RANCHC  4)  5{CTCP+BRANCH(  5 )"STCP+BRANCH(6 )  "S I NP+XTRANC 
Y1(K)=-BRANCH(1)5!SINT  +  BRANCH(2)5:COST  +  YTRAN 
Y2(K)=-BRANCH(4)5{SINT  +  BRANCH(5)J{COST  +  YTRAN 
Z1CK)=-BRANCH(1)!:CTSP-BRANCH(2);:STSP+BRANCHC3)JCC0SP-XTRANS 
Z2CK)=-BRANCH(4)j:CTSP-BRANCHC5)”STSP+BRANCHC6)5:COSP-XTRANS 
AVDIAMCK)=DIAM 
20  IF(I.EQ.NBRAN)GO  TO  30 
READC 15)  BRANCH, DI AM 
30  CONTINUE 

ENTER  TRAJECTORY  LOOP  TO  DETERMINE  MUNITION/BRANCH  ENCOUNTERS 

DO  100  1  =  1,100 
I  F(XTRAJ(  0-999999.)  135,100 

C  DETERMINE  IF  TRAJECTORY  E.XITS  SIDE  QR  BOTTOM  OF  UYER  AND 
C  DETERMINE  COORDINATES  FOR  REENTRY  PQINTS 
C  • 

135  XPLUS=BMT /SPOCP 

zout=top 

35  XOUT=XTRAJCl)  +  XPLUS 
I F(XOUT-XL IMI T )  36,36,136 

136  XOUT=0 . 

ZOUT=ZOUT-(XTRAJ(  I  )-XLIMI  T)::SPOCP 
XPLUS=CBOTTOM-ZOUT  )/SPOCP 
REPEAT= 1 
GO  TO  40 
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36  ZOUT=BOTTOM 
REPEAT=0 

ROTATE  COORDINATES  OF  TRAJECTORY  ENTRY  POINT  INTO  X'-Y’-Z' 
COORDINATE  SYSTEM 

40  YTRAJP=YTRAJ( I ) 

ZTRAJP=-XTRAJ(I  );:S  INP+ZTRAJ  ( I  )!:COSP 

ENTER  BRANCH  LOOP  AND  COMPUTE  FACTORS  USED  TO 
DETERMINE  IF  MUNITION  ENCOUNTERS  BRANCH 

KEEP=0 

DO  70  J=1,K 

AVRAD=AVDIAM(J)/2. 

FACTRF=(AVRADfFRAD):cx2 
FACTRB=CAVRAD+BRAD)5{5!2 
A=Z1(J)  -  Z2(J) 

B=Y2(J)  -  Y1(J) 

C=-Z1CJ)”B  -  Y1(J)S{A 
LENSQR=A:{}{2  +  B:{5{2 
LENGtH=SQRT(LENSQR) 

DETERMINE  IF  BRANCH  IS  PARALLEL  TO  MUNITION  TRAJECTORY  (CASE  1) 

IF(A)  53, 141 

141  I  F(B)  53,142 

142  D=(Y1CJ)-YTRAJP);!!J2+CZ1CJ)-ZTRAJP)5J”2 

IF  D  IS  GREATER  THAN  SQUARE  OF  SUM  OF  BRANCH  RADIUS  PLUS 
MUNITION  BODY  RADIUS,  SKIP  TO  NEXT  BRANCH 

IF(D-FACTRB)  143,  1*' 3, 70 

FOR  MUNITION  BODY/BRANCH  ENCOUNTER,  DETERMINE  IF  BRANCH  IS 
LARGE  ENOUGH  TO  DETONATE  MUNITION 

143  IFCAVDIAMGO-BBDIAM)  50,44,44 

DETERMINE  COORDINATES  OF  DETONATION 

44  BRDIAM=AVDIAM(J) 

YDET=YTRAJP 

ZDET=ZTRAJP 

IF(X1(J)-X2(J))  145,145,45 
145  XDET=X1(J) 


Figure  11.  (Continued) 
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GO  TO  65 

45  XDET=X2(J) 

GO  TO  65 

DETERMINE  IF  MUNITION  FUZE/BRANCH  ENCOUNTER  OCCURS,  AND  (IF 
SO)  DETERMINE  COORDINATES  OF  DETONATION 

50  I F(D-FACTRF)  44,44,70 

IF  D  IS  GREATER  THAN  SQUARE  OF  SUM  OF  BRANCH  RAUTUS  PLUS 
MUNITION  BODY  RADIUS,  SKIP  TO  NEXT  BRANCH 

53  D=((A5CYTRAJP+B5{ZTRAJP+C)/LENGTH)!{!C2 
IF(D-FACTRB)  153,153,70 

DETERMINE  IF  NORMAL  FROM  TRAJECTORY  TO  BRANCH  CENTERLINE 
INTERSECTS  BRANCH  BETWEEN  ENDS  (CASE  2) 

153  R(  1)=(Y  1(J)-YTRAJP)5!5!2+(Z  1( J)-ZTRAJP)SC!{2 
R(2)=(Y2(J)-YTRAJP)5!5{2+(Z2(J)-ZTRAJP)5C}{2 
SSQ=AMAX1(R(1),R(2))-D 
IF(SSQ-LENSQR)  154,57,57 

FOR  MUNITION  BODY/BRANCH  ENCOUNTER,  DETERMINE  IF  BRANCH  IS 
LARGE  ENOUGH  TO  DETONATE  MUNITION 

154  IF(AVDIAMCJ)-BBDIAM)  56,46,46 

DE  'ilNE  COORDINATES  OF  DETONATION 

46  SS2=R(2)-D 
S2=0. 

IFCSS2)  47,47,54 

54  S2=SQRT(SS2) 

47  S1=LENGTH  -  S2 

PX=(S2:cXl(J)  +  S1:!X2(J))/(S1  +  S2) 

55  BRDIAM=AVDIAM(J) 

YDET=YTRAJP 

ZDET=ZTRAJP 

XDET=PX 

GO  TO  65 

DETERMINE  IF  MLNITION  FUZE/BRANCH  ENCOUNTER  OCCURS,  AND  (IF 
SO)  DETERMINE  COORDINATES  OF  DETONATION 


Figure  11.  (Continued) 


38 


oooooo  ooono  oooo  r>  o  o  o  r> 


*.l.  If(l>  rACTRF)  46,4(>,70 


(  NORMAL  FROM  TRAJECTORY  TO  BRANCH  CENTERLINE  DOES  NOT  INTERSECT 

C  BRANCH  BETWEEN  ITS  ENDS  (CASE  3)  -  DETERMINE  IF  MUNITION  BODY 

C  ENCOUNTERS  ONE  OF  BRANCH  ENDS 

C 

57  MIN=1 

IF(R(2)-R(1))  164, 155, 155 
164  MIN=2 

155  IF(RCMIN)-FACTRB)  156,  156,  70 

FOR  MUNITION  BODY/BRANCH  ENCOUNTER,  DETERMINE  IF  BRANCH  IS 
LARGE  ENOUGH  TO  DETONATE  MUNITION,  AND  (IF  SO)  DETERMINE 
COORDINATES  OF  DETONATION 

156  IF(AVDIAMCJ)-BBDIAM)  63,159,  159 
159  GO  TO(59,6l),MIN 

59  PX=X1(J) 

GO  TO  55 
61  PX=X2(J) 

GO  TO  55 

DETERMINE  IF  MUNITION  FUZE/BRANCH  ENCOUNTER  OCCURS,  AND  (IF 
SO)  DETERMINE  COORDINATES  OF  DETONATION 

63  IF(RCMIN)-FACTRF)  159,159,70 

DETERMINE  Z  COORDINATES  FOR  DETONATIONS  AND  SAVE  Z  COORDINATES 
AND  BRANCH  DIAMETERS  FOR  ALL  DETONATIONS  UNTIL  ALL  BRANCHES 
IN  LAYER  HAVE  BEEN  PROCESSED 

65  ZHGT  = . 003281!{(XDET}{SINP  +  ZDETi:COSP) 

I F(ZHGT. LT. 0 .  )ZHGT=0. 

KEEP=KEEP  +  1 
KEPDET(KEEP)=ZHGT 
KEPDIAM(KEEP)=BRDIAM 
70  CONTINUE 

SET  BURST  HEIGHT  AT  MAXIMUM  Z  COORDINATE  OF  DETONATIONS,  INCRE¬ 
MENT  AIR  BURST  COUNTER,  PRINT  BURST  HEIGHT  AND  BRANCH  DIAMETER 
FOR  DETONATION  (IF  REQUIRED),  AND  FLAG  TRAJECTORY  SO  THAT  IT 
WILL  NOT  BE  CONSIDERED  IN  SUBSEQUENT  LAYERS 

IF(KEEP)  170,80 


Figure  11.  (Continued) 
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170  THEDET=KEPDET(1) 

THEDIAM=KEPDIAM(1) 

DO  71  LL=1,KEEP 

IF(KEPDETCLL)-THEDET)  71,71,171 

171  THEDET=K£PDET(LL) 

THEDIAM=KEPDIAM(LL) 

71  CONTINUE 

75  KTR  =  KTR  +  1 

I F( I  PRINT. EQ. 0)  WR1TEC6,  1004)  THEDET, THEDIAM 
XTRAJC  0=999999. 

76  DET(  I  )=THEDET 
GO  TO  100 

C 

C  SET  BURST  HEIGHT  FOR  TRAJECTORY  TO  ZERO  IF  BOTTOM  LAYER  HAS  BEEN 
C  PROCESSED  AND  IF  TRAJECTORY  DOES  NOT  HAVE  TO  BE  REENTERED 

80  I  F( BOTTOM)  85,180 

180  I  F( REPEAT)  85,181 

181  THEDET=0 . 

THEDIAM=0 . 

XTRAJC  0=999999. 

GO  TO  76 

SET  NEW  ENTRY  OR  REENTRY  COORDINATES  FOR  TRAJECTORY,  DEFINE  TOP 
OF  NEXT  LAYER,  AND  TRANSFER  TO  READ  BRANCH  DATA  FOR  NEXT  LAYER 

85  XTRAJC 0=XOUT 
ZTRAj(0=zoirr 
I  F( REPEAT)  99,100 
99  YTRAJ(I)=RANF(TRASH>'!YLIMIT 
GO  TO  35 
100  CONTINUE 
TOP=BOTTOM 
GO  TO  15 

AFTER  ALL  LAYERS  HAVE  BEEN  PROCESSED  FOR  100  TRAJECTORIES,  CALL 
SYSTEM  SUBROUTINE  BDS  TO  COMPUTE  MEAN,  VARIANCE,  AND  STANDARD 
DEVIATION  FOR  BURST  HEIGHTS,  AND  PRINT  NUMBER  OF  GROUND  BURSTS 
AND  STATISTICAL  DATA  IF  CALLED  FOR  BY  PRINT  OPTION 

110  DO  200  1  =  1,  100 
IF(DETCO)  201,200 

200  CONTINUE 

STAT(1)=STAT(5)=STAT(6)=0. 

GO  TO  210 

201  CALL.  UDSCDET,  100,STAT) 


Figure  11,  (Continued) 
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210  KTGRD  =  100  -  KTR 

lF(IPRINT.EQ.l)  GO  TO  215 

WRITE  (6,1007)  KTGRD 

WR1TE(6,1006)STAT(1),STAT(5),STAT(6) 

215  DLMMY=DUMMY+1. 

CALL  SUBROUTINE  CFRQ  TO  DETERMINE  CUMULATIVE  BURST  HEIGHT 
DISTRIBUTION  AND  INCREMENT  AZIMUTH  ANGLE.  IF  LAST  AZIMUTH  ANGLE 
HAS  NOT  BEEN  PROCESSED,  TRANSFER  TO  PROCESS  TRAJECTORIES  FOR 
NEXT  AZIMUTH  ANGLE.  OTHERWISE,  TRANSFER  TO  READ  IN  DATA  FOR 
NEXT  CASE  OR  TERMINATE  PROGRAM. 

CALL  CFRQ 
THETA=THETA+90. 

I F( THETA- 360.)  6,8,6 

1000  FORMATC 6F1 0.0, 21 2/4F10.0,  10X, I5,A1Q) 

1001  FORMmi(I4,I6,7F7.0,F6.0) 

1002  FORMAT(6F7.0,F6.0) 

1004  FORMAT(50X, 2F1 1 . 2) 

1006  FORMATC 1H0, 50X, 2 0HME AN  BURST  HEIGHT  =  ,F10.2,6H  FEET/ 

5!  6  OX,  11H  VARIANCE  =  ,F12.4/ 

::  50X, 21HSTANDARD  DEVIATION  =  ,F10.2,6H  FEET) 

1007  FORMATC  1H0, 50X,  15HTHERE  WERE  ALS0,I4,15H  GROUND  BURSTS.) 

2000  FORMATC  1H1,  60X,  4HSITE,  A10/ 1H0,  56X,  1 7HMUNITION  NUMBER  =,  1 5/1H0, 20X, 
124HMUNI TION  BODY  DIAMETER  =,F10.2,2X,  11HMILLIMETERS,  10X,  17HELEVATI 
2 ON  ANGLE  =,  F10 . 0, 2X,  7HDEGREES/1H0, 20X,  24HMUNITION  FUZE  DIAMETER  =, 
3F10 . 2, 2X,  1 1HMI LLIMETERS,  10X,  17HAZIMUTH  ANGLE  =,  F10 . 0,  2X,  7HDEGREE 

4S/1H0,  77X,  17HRANDOM  NO.  SEED  =,F12.2/1H0,20X,64FMINIMUM  BRANCH  DIA 
5METER  REQUIRED  FOR  DETONATION  GIVEN  BODY  HIT  =,F10 .2, 2X,  11FMILLIME 
6TERS/1H0, 20X,64MINIMUM  BRANCH  DIAMETER  REQUIRED  FOR  DETONATION  GIV 
.  7EN  FUZE  HIT  =,  F10.2,  2X,  1 1HMI LLIMETERS//) 

20 10  FORMATC  1H0,  54X,  5HBURST,  6X,  6HBRANCH/55X,  6HHE IGHT,  4X,  8HDIAMETER/ 

"  55X,  5HCFT . ),  6X,  5H(MM.  )/54X,  8C 1H8),  2X,  1 OC  1H5:)/) 

300  CALL  EXIT 
END 


FIGURE  11.  (Continued) 
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SUBROUTINE  CFRQ 

THIS  SUBROUTINE  DETERMINES  CUMULATIVE  BURST  HEIGHT  DISTRIBUTIONS 

DIMENSION  FREQC 150),  CFRI.QO50),  DET(IOO) 

COMMON  NMUN,  SITE,  PHI  SAV,  THETA,  X INT,  I  PUNCH,  I  PR  1  NT,  DET 

IF  THIS  IS  FIRST  OF  FOUR  CALLS  TO  THIS  SUBROUTINE,  INITIALIZE 
VARIABLES  AND  ARRAYS 

I  FCTHETA.NE  .0.0)  GO  TO  200 
DO  50  1=1,150 
50  FREQCl )=0.0 
XNUM=400 . 

DETERMINE  INTERVAL  ASSOCIATED  WITH  EACH  BURST  HEIGHT  AND  INCREMENT 
APPROPRIATE  ARRAY  LOCATION 

200  DO  208  1  =  1,  100 

D1=AINT(  100 .  *:(DET(  I )+. 005))/100. 

INC=(D1+XINT-.001)/XINT+1 .0 
FREQC !NC)=FREQ(INC)+1 . 0 
208  CONTINUE 

I  FCTHETA.NE. 270.)  RETURN 

C  IF  LAST  AZIMUTH  ANGLE  HAS  BEEN  PROCESSED,  COMPUTE,  PRINT,  AND  PUNCH 
C  (IF  REQUIRED)  HEADER  INFORMATION,  LIMITS  OF  BURST  HEIGHT  INTERVALS, 
C  AND  CUMULATIVE  PERCENTAGES  OF  DETONATIONS  IN  INTERVALS 
C 

265  CFREQC  1)=FREQ(  1)/XNUM5!100 .0 
XX=.011 
XDATA=0 . 0 
IX  =  FREQCl) 

WRI TE(6,  324) 

IF(IPRINT.EQ.O)  WRITEC6, 325) 

WRI TEC6, 261)  IX,CFREQ(1) 

I  F(  IPUNCH. EQ. 0)  GO  TO  300 
WRITEC 10, 801)  SITE, NMUN, PHISAV 
WRITEC 10,  802)  IX, CFREQC 1) 

300  DO  280  1=2,150 

CFREQC I )=CFREQ( I-1)+FREQ(  I)/XNUM”100. 

XDATA=XDATA+X I NT 
IX  =  FREQCl) 

WRITE  C 6 , 263)  XX, XDAT A,  IX, CFREQC I) 

I F( U'UNCH.EQ. 0)  GO  TO  310 

WRI  T!  (10,803)  XX,XDATA,  IX, CFREQC I) 
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:kmkiqu).i;l.9‘J.99)  go  to  1000 
280  CONTINUE 
1000  RETURN 

261  FORMAT  (1H0,  47X,  3 7H INTERVAL  NO.  OF  CUMULATIVE,  /  46X, 

"39HFROM  TO  DETONATIONS  PERCENTAGE,  /  45X, 

!{2X,  2X,  2X,  10H!!>«{!!}»{5t5«{5{,  //  45X, 

!t  13HG ROUND  BURSTS,  1 11,  FI 3. 2) 

263  FORMAT(44X,F6.2,F8.2,Ill,F13.2) 

324  FORMAT(lHO) 

325  FORMAT(lHl) 

801  FORMAT(4HSITE, A10, 5X,  15HMUNITION  NUMBER,  1 5, 5X, 15HELEVATION  ANGLE, 
1F7.2,5H  DEG.) 

802  FORMAT(20X,  1 10, FI 0.2) 

803  FORMAT(2F10. 2, 1 10, FI 0. 2) 

END 


Figure  11.  (Concluded) 
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SECTION  V 
SIMULATION  MODEL 


Tiu'  mil)  Program  uses  as  input  from  tape  the  source  and  terminal 
\,  'i  .  and  /,  coordinates  and  the'  average  diameters  of  branches  surveyed 
at  actual  forested  site's.  It  then  eiete'rmines  the  height  (or  Z  coordinate) 
lor  the’  first  branch  which  is  large  enough  to  detonate  the  munition  and 
which  is  encountered  by  randomly  selected  straight-line  trajectories 
that  pass  through  the  environment.  The  program  determines  individual 
burst  heights  for  100  munition  trajectories  from  each  of  four  azimuth 
angles  (i.e.,  0,  90,  180,  and  2 70  degrees)  and  for  any  elevation  angle 
specified  as  input  and  then  prepares  a  cumulative  burst  height  distribution 
table  for  the  400  trajectories. 


TABLES  AND  ARRAYS 

No  taiiles  are  required  in  the  PHD  Program,  and  the  arrays  are 
used  only  for  storage  of  computed  data.  The  DIMENSION  and  COMMON 
statements  used  in  Program  MAIN  are: 


DIMENSION  XTRAJ(100),YTRAJ(  100), ZTRAJ(IOO),  BRANCH(6),  R(2),  Xl(  1500) 
1,  Y1(1500),Z1(  1500),  X2(  1500),  Y2(  1500), Z2C1500), AVDIAMC1500),  DETC100 
2 )  ,  KEPDET( 3 0 0  ) ,  KEPD I  AM(  3 0 0  ) 

DIMENSION  STATC8) 

COMMON  NMUN,  SITE, PHI  SAV,  THETA, XINT,  I  PUNCH,  IPRINT,DET 
COMMON  XI, Yl,  Z1,X2,  Y2,  Z2, AVDIAM 


and  the  DIMENSION  and  COMMON  statements  used  in  Subroutine  CERQ 
a  re: 


DIMENSION  FREQC150),  CFREQCI30),  DET(100) 

COMMON  NMUN,  S I TE  ,  PH  I S  AV,  THETA,  X I  NT,  I  PUNCH,  I  PR  I  NT,  DET 


The  COMMON  variables  used  in  Program  MAIN  and  not  in 
Subroutine  CERQ  are  placed  only  in  Program  MAIN  to  reduce  core 
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storage  requirements.  The  system  loader  entry  point  is  placed  at  the 
beginning  of  BLANK  COMMON.  BLANK  COMMON  is  larger  than  the 
loader,  and  the  core  storage  requirements  are  therefore  reduced  by  the 
size  of  the  loader. 

All  variables  which  appear  in  these  DIMENSION  and  COMMON 
statements  are  individually  defined  in  the  List  of  Symbols  and  Abbre¬ 
viations  (Simulation  Model). 


DISCUSSION  OF  SIMULATION 


Program  MAIN 


Program  MAIN  computes  the  munition  burst  heights  for  each  of 
the  munition  trajectories  being  considered.  The  first  executable 
statements: 


CALL  EOF(KEOF) 

8  ASSIGN  300  TO  KEOF 

READC5,  1000)  PHI, XLIMIT,YLIMIT,ZLIMIT, DUMMY, XINT,  I  PUNCH,  I  PRINT, 
1 BD I  AM,  FDI  AM,  BBD I  AM,  BFD I  AM,  NMUN  ,  S I TE 


provide  a  transfer  to  terminate  the  program  if  an  end-of-file  is  read 
on  the  input  file,  and  they  read  the  two  data  cards  defining  parameters 
for  the  case.  EOF  is  a  system  routine  which  allows  the  program  to 
retain  control  after  an  end-of-file  has  been  read.  The  next  statements: 


PHISAVE=PHI 

THETA=0. 

6  REWIND  15 
PH I =PH I  SAVE 
DO  1  1=1,100 

1  DET(I)=0. 

DO  2  1=1,300 
KEPDIAMC I )=0 . 

2  KEPDET(I)=0. 
ASSIGN  110  TO  KEOF 
KEEP=0 
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initialize  several  program  variables  and  arrays,  rewind  the  binary  tape 
containing  the  branch  input  data,  and  provide  a  transfer  to  perform  the 
statistical  calculations  if  an  end-of-file  is  read  op  the  branch  input  data 
tape.  The  next  statements: 


I  FCTHETA.NE.  0.  .ANDIPRINT.EQ.  1)  GO  TO  7 

WRITE(6, 2000)SITE,NMUN,BDIAM,  PHI,  FDI  AM,  THETA,  DUMMY,  BBDI  AM, 
1BFDIAM 

IF(IPRINT.EQ.l)  GO  TO  7 
WRITE  (6,2010) 


write  the  header  information  for  the  case  on  the  output  file  as  called  for 
by  the  print  option.  The  next  statements: 


7  KTR^O 
PHISAV  =  PHI 
PHI=-PHI 
XTRAN=0. 

YTRAN=0. 

I THETA=THETA/ 90 .  +  1. 

GO  TO(5, 3, 3, 4), ITHETA 

3  IFCITHETA.EQ.3)  XTRAN=XLIMIT 
YTRAN=YLIM1T 

GO  TO  5 

4  XTRAN=XLIM1T 

5  FRAD=FDIAM/2. 

BRAD=BDIAM/2. 

AZ I  MLfTH=  .01745  3295"THETA 
PH  I  = .  0 1  74532955!PH1 
SINT=SIN(AZIMUTH) 
SINP=S1N(PHI) 

COST=COS (AZIMUTH) 
COSP=COS(PHI) 

CTCP=COST::COSP 

CTSP=C0ST”S1NP 

STSP=SINT;CSINP 

STCP=SINT"COSP 

XTRAMC=XTRAM”COSP 

XTRANS=XTRAN”SIWP 

SPOCP^S 1 NP/COSP 


perform  initial  calculations  which  are  used  later  in  the  program.  These 
calculations  include  changing  the  sign  of  the  elevation  angle,  determining 
the  translation  values  required  to  place  all  transformed  branch  coordi¬ 
nates  into  the  first  octant,  computing  the  munition  fuze  and  body  radii, 
and  computing  trignometric  functions  of  the  azimuth  and  elevation  angles. 
The  next  statements: 


XLIMO=XLIMI T/10 . 
YLIMO=YLIMI T/10 . 
XLIMOO=XLIMIT/100. 
YLIM00=YLIMI T/100 . 

CALL  RANSETCDimO 
DO  10  1  =  1,  10 
XCELL=XLIMO::(  I  -1 ) 
K=(I~1)::10 
DO  10  d=l,  10 
YCELL=YLIM0”(J-1) 

NRNXY= RAN  FC  T  RASH  )5: 1 0  0 . 

NRNX=NRNXY/10 

RNX=NRNX 

RNY=NRNXY-(NRNX!{10) 

L=K+J 

XTRAJ(L)=XC£LL+RNX"XLIM00 
YTRAJ(L)=YCELL+RNY"YLIM00 
10  ZTRAJ(L)=ZLIMIT 


initialize  the  system  uniform  random  number  generator,  subdivide  the 
X-Y  plane  into  100  rectangles,  and  randomly  select  munition  trajectory 
X  and  Y  entry  coordinates  which  pass  through  each  of  the  rectangles 
(Equations  1,  2,  and  3).  The  sequential  number  of  the  layer  of  branch 
input  data,  the  number  of  branches  in  the  layer,  the  Z  coordinate  at  the 
bottom  of  the  layer,  and  the  X,  Y,  and  Z  coordinates  and  average  diam¬ 
eter  of  the  first  branch  are  read  in  by  the  statements: 


TOP=ZLIMIT 

15  READC15)  LEVEL,  NBRAN,  BOTTOM,  BRANCH,  DI  AM 
K=0 

BMT =B0TT0M-T0P 
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and  the  IX)  LOOP: 


DO  30  I  =  1,NBRAN 
I F(DIAM.LT. BFDIAM)GO  TO  20 
K=K  +  1 

X1(K)=BRANCH(1)!{CTCP+BRANCHC2)“STCP+BRANCH(3):cSINP+XTRANC 
X2(K)=BRANCH(4)5{CTCP+BRANCH(5)!!STCP+BRANCHC6)!{SINP+XTRANC 
Y1(K)=-BRANCH(1)"SINT  +  BRANCH(2)S!COST  +  YTRAN 
Y2(K)=-BRANCH(4)!!SINT  +  BRAN  CH(  5)  “COST  +  YTRAN 
Zl(K)=-BRANCH(l)54:TSP-BRAN(>l(2)s!STSP+BRANCH(3)“COSP-XTRANS 
Z2(K)=-BRANCH(4)“CTSP-BRANCH(5)!{STSP+BRANCH(6)“COSP-XTRANS 
AVDIAM(K)=DIAM 
20  IFCI  .EQ.NBRAN)GO  TO  30 
RE  ADC  15)  BRANCH, DIAM 
30  CONTINUE 


l  eads  in  the  source  and  terminal  X,  Y,  and  Z  coordinates  of  the  re¬ 
maining  branches,  eliminates  those  branches  with  a  diameter  less  than 
the  minimum  branch  diameter  that  will  detonate  the  munition  if  en¬ 
countered  by  the  munition  fuze,  and  transforms  the  X,  Y,  and  Z  co¬ 
ordinates  of  the  branch  source  and  terminal  ends  into  a  new  X'-Y'-Z’ 
coordinate  system  with  the  X'  axis  parallel  to  the  paths  bf  the  munition 
trajectories  (Equations  4,  5,  and  6).  The  branch  coordinates  are  rotated 
through  the  azimuth  angle  about  the  Z  axis,  translated  into  the  first 
octant,  and  rotated  through  the  elevation  angle  about  the  new  Y*  axis. 

The  trajectory  loop  is  entered  to  determine  the  munition/ branch 
encounters,  and  the  statements: 


DO  100  1=1,100 
I FCXTRAJC  D-999999.)  135,  100 


cause  the  program  to  proceed  to  the  next  munition  trajectory  if  a  flag 
indicates  that  a  detonation  has  been  identified  for  the  munition  trajectory 
m  a  previous  layer.  The  next  statements: 


4  H 


135  XPLUS=BMT /SPOCP 
ZOUT=TOP 

35  XOUT=XTRAJ(  I  )  +  XPLUS 
IFCXOITT-XUMIT)  36,36,136 


compute  the  X  coordinate  of  the  point  at  which  the  munition  trajectory 
intersects  the  plane  defining  the  bottom  of  the  layer  (Equation  7).  If 
the  munition  trajectory  exits  the  layer  on  the  side  rather  than  on  the 
bottom,  the  reentry  coordinates  for  the  munition  trajectory  are  com¬ 
puted  by  the  statements  (Equations  9  and  10): 


136  XOUT=0. 

ZOUT=ZOlfT-(XTRAJ(  I  )-XLIMI  T)5!SP0CP 
XPLUS=CB0TT0M-Z0UT)/ SPOCP 
REPEAT  = 1 
GO  TO  40 
36  ZOlfT=BOTTOM 
REPEAT=0 


If  the  munition  trajectory  exits  through  the  bottom  of  the  layer,  the  re¬ 
entry  coordinates  for  the  next  layer  are  computed  (Equations  13,  14, 
and  15).  The  reentry  coordinates  are  then  transformed  so  that  the 
munition  trajectory  paths  ar^  parallel  to  the  X'  axis  (Equations  17  and 
18)  by  the  statements: 


40  YTRAJP=YTRAJ(I) 

ZTRAJP=-XTRA*J(  I  )”S INP+ZTRAJC I  )!:C0SP 


The  branch  loop  is  then  entered,  and  the  statements: 


KEEP=0 

DO  70  J=1,K 

AVRAD=AVDIAM(J)/2. 

FACTRF=(AVRAD+FRAD),{}{2 

FACTR8=(AVRAD+BRAD)S{5{2 
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A=Z1(J)  -  72CJ) 
B-Y2CJ)  -  Y1(J) 
C=-Z1(J)::B  -  Y1(J)::A 
LENS(}R=A::::2  +  Br:2 

length=sqrtClen^.qr) 


compute  several  factors  which  are  used  later  in  the  propram  to  determine 
if  the  munition  body  or  fuze  encounters  a  branch  (Equations  33  through 
36).  The  next  statements: 


I F(A)  53,141 
141  I F(B)  53,142 


check  to  determine  if  both  A  and  B  are  zero  (i.  e.  ,  the  branch  centerline 
is  parallel  to  the  munition  trajectory).  If  so,  the  requirements  for  Case  1 
are  satisfied,  and  the  square  of  the  distance  in  the  transformed  Y'-Z' 
plane  between  the  munition  trajectory  and  the  branch  centerline  (D)  is 
computed  by  the  statement  (Equation  22): 


142  D=(Y1(J)-YTRAJP)”"2+(Z1(J)-ZTRAJP):":2 


The  value  of  D  is  then  compared  to  the  square  of  the  sum  of  the  munition 
body  radius  and  the  branch  radius  (FACTRB)  by  the  statement  (Inequality 


23): 


IF(D-FACTRB)  143,143,70 


If  D  is  larger  than  FACTRB,  the  munition  body  does  not  encounter  the 
branch,  and  the  program  proceeds  to  the  next  trajectory/ branch  com¬ 
bination.  If  D  is  smaller  than  FACTRB,  the  munition  body  encounters 
the  branch,  and  the  statement  (Inequality  25): 


143  I  F(AVDIAM(J)-BBDIAM)  50,44,44 


determines  if  the  branch  is  large  enough  for-  detonation  by  the  munition 
body.  If  the  branch  diameter  is  larger  than  the  minimum  diameter 
required  for  detonation  by  the  munition  body,  the  statements  (Equations 
26,  27,  and  28): 


44  BRDIAM=AVDIAM(J) 
YDET=YTRAJP 
ZDET=ZTRAJP 

IF(X1(J)-X2(J))  145,145,45 
145  XDET=X1(J) 

GO  TO  65 

45  XDET=X2(J) 

GO  TO  65 


compute  the  transformed  X',  Y\  and  7*  coordinates  of  the  munition 
detonation.  If  the  branch  diameter  is  smaller  than  the  minimum  diam¬ 
eter  required  for  detonation  by  the  munition  body,  the  statement 
(Inequality  24): 


50  JF(D-FACTRF)  44,44,70 


determines  whether  the  munition  fuze  encounters  the  branch.  If  the 
munition  fuze  encounters  the  branch,  the  transformed  X',  Y\  and  Z* 
coordinates  for  the  detonation  are  computed  as  described  earlier.  If 
the  munition  fuze  does  not  encounter  the  branch,  the  program  proceeds 
to  the  next  trajectory/branch  combination. 

If  the  values  of  A  and  B  are  not  zero,  the  branch  centerline  is  not 
parallel  to  the  munition  trajectory,  and  the  statements: 


53  D=((A55YTRAJP+B!{ZTRAJP+C)/LENGTH)5C}S2 
IF(D-FACTRB)  153,  153.70 


compute  the  square  of  the  normal  in  the  transformed  Y ‘ -7l  plane  from 
the  tmimtion  trajectory  to  the  branch  centerline  (Equation  35)  and  com¬ 
pare  tins  value  (D)  to  the  square  of  1  he  sum  of  the  munition  body  radius 
and  ttie  branch  radius  (FACTUM).  If  U  is  larger  than  FACTRB  (Inequality 
L’3),  the  munition  body  does  not  encounter  the  branch,  and  the  program 
proceeds  to  the  next  trajectory /branch  combination.  If  D  is  smaller  than 
FACTRB,  the  statements  (Equations  31  and  32  and  Inequality  30): 


153  R(  l)=(Yl(J)-YTRAuP)!{::2+(Zl(J)-ZTRAJP)::"2 
R(2)=(Y2(j)-YTRAJP):::i2+(Z2(J)-ZTRAJP)::5!2 
SSQ=AMAX1(R(  0/  R(  2))-D 
IF(SSQ-LENSQR)  154, 57, 57 


determine  whether  the  normal  in  the  transformed  Y'-Z'  plane  from  the 
munition  tra)ectory  to  the  branch  centerline  intersects  the  branch  center- 
line  between  the  branch  ends.  If  so,  the  requirements  for  Case  2  are 
satisfied,  and  the  statement  (inequality  25): 


154  IFCAVDIAM(J)-BBDIAM)  56,46,46 


de  ermines  if  the  branch  is  large  enough  for  detonation  by  the  munition 
body.  If  the  branch  diameter  is  larger  than  the  minimum  diameter 
required  for  detonation  by  the  munition  body,  the  statements  (Equations 
37  through  41): 


46  SS2=R(2)-D 
S2=0 . 

1FCSS2)  47,47,54 

54  S2=SQRT(SS2) 

47  S1=LENGTH  -  52 

PX=(S2;{X1(J)  +  S1”X2CJ))/(S1  +  52) 

55  BRDIAM=AVDIAM(j) 

YDET=YTRAJP 

ZDET=ZTRAJP 

XDET-PX 

GO  TO  65 


.  I 


compute  the  transformed  X' ,  Y* ,  and  z'  coordinates  of  the  munition 
detonation.  If  the  branch  diameter  is  smaller  than  the  minimum  diam¬ 
eter  required  for  detonation  by  the  munition  body,  the  statement 
(Inequality  24): 


56  IF(D-FACTRF)  46,46,70 


deteri  n.  es  whether  the  munition  fuze  encounters  the  branch.  If  the 
munition  fuze  encounters  the  branch,  the  transformed  X',  Y' ,  and  Z' 
coordinates  are  computed  as  described  earlier.  If  the  munition  fuze 
does  not  encounter  the  branch,  the  program  proceeds  to  the  next 
trajectory  /  branch  combination. 

If  the  normal  in  the  transformed  Y' -Z!  plane  from  the  munition 
trajectory  to  the  branch  centerline  does  iot  intersect  the  branch  center¬ 
line  between  the  branch  ends,  the  requi,  merits  for  Case  3  are  satisfied, 
and  the  statements  (Inequality  43): 


57  MIN=1 

IF(R(2)-R(1))  164,  155,155 
164  MIN=2 

155  IF(RCMIN)-FACTRB)  156,156,70 


determine  if  the  munition  body  encounters  either  of  the  branch  ends. 
If  the  munition  body  does  not  encounter  either  of  the  branch  ends,  the 
program  proceeds  to  the  next  trajectory/branch  combination.  If  the 
munition  body  encounters  a  branch  end,  the  statement  (Inequality  25): 


156  IF(AVDIAMCd)-BBDIAM)  63,159,159 


determines  if  the  branch  is  large  enough  for  detonation  by  the  munition 
bod> .  If  the  branch  diameter  is  larger  than  the  minimum  diameter 
required  for  detonation  by  the  munition  body,  the  statements: 


159  GO  T0(59,61),MIN 
59  PX=X1(J) 

GO  TO  55 
61  PX=X2(J) 

GO  TO  55 


select  the  transformed  X'  coordinate*  and  compute  the  Y*  and  /'  co¬ 
ordinates  of  the  munition  detonation  (Equations  46,  46,  and  47).  If  the 
branch  diameter  is  smaller  than  the  minimum  diameter  required  for 
detonation  by  the  munition  body,  the  statement  (Inequality  44): 


63  IF(RCMIN)-FACTRF)  159,159,70 


determines  whether  the  munition  fuze  encounters  the  branch.  If  the 
munition  fuze  encounters  the  branch,  the  transformed  X',  Y',  and  Z' 
coordinates  for  the  detonation  are  computed  as  described  earlier.  If 
the  munition  fuze  does  not  encounter  the  branch,  the  program  proceeds 
to  the  next  trajectory / branch  combination. 

In  the  last  statements  of  the  branch  loop: 


65  ZHGT=. 00 32 81::(XDET::SINP  +  ZDET::COSP) 
I  FCZHGT . LT. 0 . )ZHGT-0. 

KEEP=KEEP  +  1 
KEPDET(KEEP)-ZHGT 
KEPDIAM(KEEP)=BRDIAM 
70  CONTINUE 


the  Z  coordinate  for  the  detonation  is  computed  (Equation  29),  and 
the  Z  coordinates  and  corresponding  branch  diameters  are  saved  for 
all  detonations  until  all  branches  in  the  layer  have  been  processed. 

If  any  detonations  occur  in  the  layer,  the  statements: 


IF(KEEP)  170/ SO 

170  THEDET=KEPDETU) 
THEDI/W=KEPDIAM(1) 

DO  "1  LL=1/KEEP 

I  F(KEPDET(LL)-THEDET)  71/71/171 

171  THEDET=KEPDET(LL) 
THEDIAM=KEPDIM(LL) 

71  CONTINUE 


set  the  burst  height  for  the  trajectory  at  the  maximum  Z  coordinate  of 
the  detonations.  The  next  statements: 


75  KTR  =  KTR  +  1 

IF(IPRINT.EQ.O)  WRITEC6/1004)  THEDET,THEDIAM 
XTRAJC  0=999999. 

76  DET( I )=THEDET 
GO  TO  100 


increment  the  counter  for  the  number  of  air  bursts,  print  the  burst 
height  and  the  branch  diameter  which  caused  the  detonation  if  called 
for  by  the  print  option,  and  flag  the  munition  trajectory  so  that  it  will 
not  be  considered  in  the  subsequent  layers  of  branches.  If  no  detonations 
occur  for  the  munition  trajectory  in  the  layer  of  branches,  the  state¬ 
ments  : 


80  I F( BOTTOM)  85,180 

180  I F(REPEAT)  85,181 

181  THEDET=0. 
THEDi^M=0. 

XTRAJC  0=999999. 
GO  TO  76 


set  the  burst  height  for  the  munition  trajectory  to  zero  (i.  e.  ,  a  ground 
burst)  if  the  bottom  layer  of  branches  has  just  been  processed  and  if  the 
munition  trajectory  does  not  have  to  be  reentered  into  the  layer.  Other¬ 
wise,  the  statements: 
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85  XTRAJC 1 )=XOUT 
ZTRAJC 1 )=ZOUT 
i F( REPEAT)  99,100 
99  YTRAJC I )=RANF(TRASH)"YLIMI T 
GO  TO  35 
100  CONTINUE 


set  the  entry  or  reentry  coordinates  to  the  previously  computed  values. 
These  will  be  either  the  entry  coordinates  for  the  next  layer  or,  if  the 
munition  trajectory  did  not  exit  the  layer  through  the  bottom,  the  re¬ 
entry  coordinates  for  the  current  layer.  After  the  100  munition  tra¬ 
jectories  have  been  processed  for  the  layer  of  branches,  the  statements: 


T0P= BOTTOM 
GO  TO  15 


define  the  top  of  the  next  layer  and  provide  a  transfer  to  read  the  branch 
data  for  the  next  layer.  After  all  layers  have  been  processed  for  the  100 
munition  trajectories,  the  statements: 


110  DO  200  1=1,100 
IFCDETCO)  201,200 

200  CONTINUE 

STATC 1 )=STAT(5)=STAT(6)=0 . 
GO  TO  210 

201  CALL  BDS(DET, 100,STAT) 


rail  System  Subroutine  i-iDS  to  compute  the  mean,  variance,  and 
standard  deviation  for  the  burst  heights  of  the  100  munition  trajectories 
if  any  air  bursts  have  occurred.  The  next  statements: 


5  G 


210  KTGRD  =  100  -  KTR 

IF(lPRINT.EQ.l)  GO  TO  215 
WRITE  (6,1007)  KTGRD 
WRITE(6, 1006)STAT( 1),STAT(5),STAT(6) 
215  DUMMY =DUMMY  + 1 . 


print  the  number  of  ground  bursts  and  the  statistical  parameters  if 
called  for  by  the  print  option.  Finally,  the  statements: 


CALL  CFRQ 
THETA=THETA+90. 
IFCTHETA-360.)  6,8,6 


call  Subroutine  CFRQ  to  determine  the  cumulative  distribution  of  the 
burst  heights  and  increment  the  azimuth  angle.  If  all  azimuth  angles 
have  not  been  processed,  a  transfer  is  made  to  process  100  munition 
trajectories  for  the  next  azimuth  angle.  Otherwise,  a  transfer  is  made 
to  read  in  data  for  the  next  case.  The  statements: 


300  CALL  EXIT 
END 


terminate  the  program. 


Subroutine  CFRQ 

The  purpose  of  this  subroutine  is  to  compute,  print,  and  optionally 
punch  the  cumulative  burst  height  distribution  for  the  munition  and  ele¬ 
vation  angle  being  considered.  The  first  executable  statements: 


IF(THETA.NE.O.O)  GO  TO  200 
DO  50  1=1,150 
50  FREQ( I )=0 . 0 
XNUM=400. 


initialize  the  cumulative  burst  height  distribution  array  and  set  the  total 
number  of  data  points  if  this  call  to  t he  subroutine  is  for  the  first  azi¬ 
muth  angle  (i.e.,  0  degrees).  The  DO  LOOP: 


200  DO  208  1  =  1,100 

D1=AINT(100.::(DET(I)+.005))/100. 
INC=(D1+XINT-.001)/XINT+1 .0 
FREQC INC)=FREQ( INC)+ 1 . 0 
208  CONTINUE 


determines  the  interval  associated  with  each  burst  height  and  increments 
the  appropriate  array  location.  If  the  100  munition  trajectories  for  the 
last  azimuth  angle  (i.e.,  270  degrees)  are  being  processed,  the  state¬ 
ments: 


1 FC THETA. NE. 270.)  RETURN 
265  CFREQ(1)=FREQ(1)/XNUM"100.0 
XX=.011 
XDATA=0 . 0 
IX  =  FREQC 1 ) 

WRITEC6, 324) 

I F( I  PRINT. EQ. 0)  WRITEC6, 325) 
WRITE(6,26i)  1X,CFREQ(1) 
IF(IPUNCH.EQ.O)  GO  TO  300 
WRI TE(10, 801)  SITE,NMUN,PHISAV 
WRITEClO, 802)  IX,CFREQC 1) 


print  the  header  information  and  compute  and  print  the  limits  of  the  first 
burst  height  interval,  the  number  of  detonations  in  the  interval,  and  the 
cumulative  percentage  of  detonations  in  the  interval  (Equation  48). 

These  values  are  also  punched  (or  written  on  Logical  File  10)  if  called 
for  by  the  punch  option.  Finally,  the  statements: 


300  DO  280  1=2,150 

CFREQC I  )=CFREQ(  I  - 1  )+FREQ(  I  )/XNUM::l  00 . 
XDATA=XDAT  Af X I NT 
IX  =  FREQC I) 


>8 


310 

280 

1000 


WRITE  (6,263)  XX,  XDATA,  IX,CFREQ( I ) 
IF(IPUNCH.EQ.O)  GO  TO  310 
WRITEClO, 803)  XX , XDAT A, 1 X, CFREQC I ) 
XX-XX+XINT 

I F(CFREQ( I ) .GE . *)‘) . 99)  GO  TO  1000 

CONTINUE 

RETURN 

END 


compute,  print,  and  optionally  punch  the  remaining  burst  height  inter¬ 
vals,  numbers  of  detonations,  and  cumulative  percentages  of  detonation 
and  stoj)  when  the  cumulative  burst  height  value'  reaches  100  percent. 
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